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ABSTRACT 


Intensities of EUV spectral lines were measured as a function of radius off the 
solar limb by two flights of Goddard’s Solar EUV Rocket Telescope and Spectro- 
graph (SERTS) for three quiet sun regions. The density scale height, line-ratio 
densities, line-ratio temperatures, and emission measures were determined. The 
line-ratio temperature determined from the ionization balances of Arnaud and 
Rothenflug (1985) were more self-consistent than the line-ratio temperatures ob- 
tained from the values of Arnaud and Raymond (1992). Limits on the filling 
factor were determined from the emission measure and the line-ratio densities for 
all three regions. 

The relative abundances of silicon, aluminum, and chromium to iron were de- 
termined. Results did agree with standard coronal relative elemental abundances 
for one observation, but did not agree for another. Aluminum was overabundant 
while silicon was underabundant. 

Heating was required above 1.15 solar radii for all three regions studied. For 
two regions, local nonconductive heating is needed for any filling factor, and in all 
three regions for filling factor of 0.1. 


iii 


PmOKNNt FAGS BLANK NOT FILM0P 




Contents 


1 Introduction 1 

1.1 Overview of Solar Structure 1 

1.2 Features 3 

1.3 Solar Emission 5 

1.4 Solar Corona g 

1.4.1 Heating g 

1.4.2 Elemental Abundance Variation g 

1.4.3 Magnetic Field Structures 7 

1.5 Outline of Dissertation 7 

1.6 Acknowledgment 7 

2 Theory g 

2.1 Density Scale Height and the Scale-height Temperature 9 

2.2 Coronal Elemental Composition 10 

2.3 The Atomic Physics of Coronal Emission Lines 12 

2.3.1 Collisional Excitation of Coronal Ions 14 

2.3.2 Ionization Balance 15 

2.4 The Line of Sight Intensity 15 

2.5 Spectroscopic Diagnostics lg 

2.5.1 Density-Sensitive Lines lg 

2.5.2 Electron Temperature Determination 19 

2.5.3 Emission Measure 20 

2.5.4 Relative Elemental Abundance 21 

2.6 Line Emissivity Effect on Scale Height 21 

3 Instrument 23 

3.1 The SERTS Instrument 23 

3.2 Data Reduction 24 

3.3 Observations 30 

3.3.1 SERTS 1989 ^ ’ ’ . ’ ’ ’ ’’’’’’’’’ ’ 30 

3.3.2 SERTS 1991 30 

4 Analysis 33 

4.1 Line Selection and Determination of Spectral Line Intensity 33 

4.1.1 SERTS 1991 ’ ’ 34 

v 


PAGK BLANK NOT FILMED 



4.1.2 SERTS 1989 

4.1.3 Sets of Ionization Balance Calculations 

4.2 Scale Height 

4.3 Line-Ratio Electron Density 

4.4 Line-Ratio Temperature 

4.4.1 SERTS 1991 

4.4.2 SERTS 1989 

4.4.3 Comparison of Different Ionization Balance Calculations 

4.4.4 Comparison of Line-Ratio Temperatures and Scale Height 

4.5 Emission Measure 

4.6 Relative Elemental Abundances 


35 

35 

35 

37 

37 

37 

41 

46 

46 

50 

56 


5 Implications of the Coronal Parameters 

5.1 Filling Factor and Minimum Electron Density 

5.2 Heating 

6 Upper Limits on Bi-Maxwellian Electrons 

6.1 Symbols 

6.2 Observations 

6.3 Theoretical Limit 

7 Conclusions 

7.1 Relative Elemental Abundances 

7.2 Radial Dependence of Coronal Parameters . 

7.3 Heating 

7.4 Bi-Maxwellian 


65 

65 

66 

71 

72 
72 
75 

79 

79 

79 

80 
80 


8 Future Applications of SERTS and SOHO 

8.1 SERTS 

8.2 SOHO CDS/NIS 


A Solar Plasma Parameters 


83 


List of Figures 

1.1 The solar cross section 2 

1.2 Yohkoh X-ray image of the Sun 4 

2.1 Relative Elemental Abundances of the Corona 11 

2.2 Illustration of Line of Sight 16 

2.3 Theoretical Fe XIII Line-Ratio Density 19 

3.1 The SERTS Instrument 24 

3.2 Layout of SERTS Spectrograph 25 

3.3 SERTS Spectrum 26 

3.4 Wide-Lobe SERTS Spectrum 27 

3.5 SERTS 1989 Pointing Positions 28 

3.6 SERTS 1991 Pointing Positions 29 

3.7 Sunspots Maps for SERTS 1989 31 

3.8 Sunspots Maps for SERTS 1991 32 

4.1 Edges of the Lobes 34 

4.2 Contribution Function 36 

4.3 Scale Height 38 

4.4 Line-Ratio Temperature using IBAL1 42 

4.5 Line-Ratio Temperature using IBAL2 44 

4.6 IBALl Emission Measures 52 

4.7 IBAL2 Emission Measures 54 

4.8 Relative Elemental Abundances using IBALl 59 

4.9 Relative Elemental Abundances using IBAL2 60 

4.10 Aluminum Abundance for SERTS 1989 62 

5.1 Heating Constraints 68 

6.1 Electron Distribution 72 

6.2 Ionization Cross Section 73 

6.3 Theoretical Value of the Log of the Iron Intensity Ratio vs the Log of the Temperature. 74 

6.4 Theoretical Upper Limits on a 78 


vii 



List of Tables 


4.1 Spectral Line Information 

4.2 Temperature from Scale Height 

4.3 Line-Ratio Temperature From SERTS 1991 Southern Edge . 

4.4 Line-Ratio Temperature From SERTS 1989 Southern Edge . 

4.5 Line-Ratio Temperature From SERTS 1989 Southern Edge . 

4.6 Comparison Scale Height vs. Line-Ratio Temperature . . . . 

4.7 Emission Measure 

4.8 Relative Elemental Abundances 

5.1 Minimum Electron Density and Filling Factor 

5.2 Comparison of Radiation to the Divergence of the Heat Flux 

6.1 Values of Parameter /? 

6.3 Temperature 

6.4 Mixing Fraction, a 

6.2 Values of Parameter — /? Fe ) 


36 

40 

47 

48 

49 
51 
57 
63 

67 

67 

75 

76 
76 
76 



List of Abbreviations 


Abbreviation Meaning 

CDS/NIS Coronal Diagnostic Spectrometer/Normal Incidence Spectrometer 

CDS/GIS Coronal Diagnostic Spectrometer/Grazing Incidence Spectrometer 

EUV Extreme-ultraviolet light 

FIP First Ionization Potential 

GSFC Goddard Space Flight Center 

IBAL1 First Set of Ionization Balances 

IBAL2 Second Set of Ionization Balances 

ISAS Institute of Space and Astronautical Science 

IR Infrared 

MSIS-86 Mass Spectrometer and Incoherent Scatter 
NASA National Aeronautics and Space Administration 

SEP Solar Energetic Particles 

SERTS Solar Extreme-ultraviolet Rocket Telescope and Spectrograph 

SOHO Solar and Heliospheric Observatory 

SOLEX Solar X-rays 

SURF-II Synchrotron Ultraviolet Radiation Facility 


IX 



Chapter 1 


Introduction 


Unlike other stars, the Sun is close enough for 
us to spatially resolve its features. This allows 
us to do a more detailed comparison of measure- 
ment and modeling of the Sun than is possible 
with other stars. The study of the Sun involves 
a wide range of physics, including atomic physics, 
nuclear physics, and plasma physics, all of which 
are required to fully understand solar phenom- 
ena. Since a variety of conditions exist on the 
Sun which cannot be reproduced in a terrestrial 
laboratory, the Sun offers a unique opportunity 
to test theories for the corresponding physical 
processes [45]. A variety of interesting questions 
about the Sun are now being studied, including: 
the method by which energy is transferred out- 
ward from the convection zone into the upper 
solar atmosphere, the source and nature of varia- 
tions in the solar wind, solar activity, flares, and 
variations of elemental abundances from region 
to region on the Sun. 

The work in this dissertation is based upon 
analysis of data from the Solar Extreme- 
ultraviolet Rocket Telescope and Spectrograph 
(SERTS), an instrument built by members of 
the Solar Physics Branch of National Aeronau- 
tics and Space Administration (NASA)’s God- 
dard Space Flight Center (GSFC), and flown 
on sounding rockets in 1989 and 1991 [43, 55]. 
SERTS obtains in first order Extreme-ultraviolet 
(EUV) (250 A to 450 A) images and spectra of 
different regions of the Solar corona. From im- 
aged spectra of the solar corona off the solar disk, 
the relative elemental abundance of, and heating 
constraints on, the observed corona have been 


determined. 

Relative elemental abundances in the solar 
corona have, over the last decade, been found to 
differ from those in the photosphere. Recently, 
variations in relative elemental abundance from 
region to region in the corona have been shown to 
exist. These variations are important for several 
reasons, which include their effects on coronal di- 
agnostics, and knowing them can lead to a bet- 
ter understanding of where the coronal material 
originates. 

To put these results into context, we will first 
provide some general background information 
on the Sun. This will include a description of 
the general structure of the Sun and of features 
on the solar disk, as well as description of the 
methods with which some different parts of the 
electro-magnetic spectrum are used to study the 
Sun. We will also provide a more in-depth dis- 
cussion of the solar corona. 

1.1 Overview of Solar 
Structure 

The Sun is not homogeneous, but instead con- 
tains many physically distinct regions. Each re- 
gion will be discussed separately. The boundaries 
of these regions are somewhat arbitrary, but for 
the purpose of this dissertation the exact bound- 
ary locations are not important (Figure 1.1 for 
solar cross section [11]). To fully understand any 
one region, the effects on that region by other 
regions need to be understood. The following 
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Figure LI. The solar cross section. Interior lay- 
ers of the sun, with the solar corona in the back- 
ground [11]. 

brief description of the principal solar regions is 
greatly expanded in a recent text by Stix [52]. 
See Appendix A for typical plasma parameters 
for photosphere, chromosphere, and corona and 
some useful equations. 

The innermost solar region is the central core, 
where the energy output of the Sun is produced 
by the fusion of four hydrogen nuclei into a 
helium nucleus in a chain of reactions. The 
core density at Sun center is on the order of 
100 gm cm -3 and the temperature is more than 
10 7 K* The core extends outward beyond the re- 
gion of nuclear reactions, through a zone where 
energy is transferred radiatively by photons. 

Above the core lies the convection zone, where 
the temperature has dropped sufficiently for re- 
combination of helium and hydrogen to occur, 
increasing the plasma opacity. This, in turn, in- 
creases the temperature gradient, so convection 
instead of radiative heat transfer becomes the 
primary mode of energy transport. The convec- 
tive motions produce a variety of waves which are 
transmitted upward out of this layer. Through- 
out this region, the temperature decreases with 
height from the order of several times 10 5 K to 
less than 10 4 A\ and the density decreases with 
height from 0.1 to 0.01 gm cm~ 3 . The solar mag- 
netic field is believed to be generated at the base 
of this region, and then to rise and emerge from 


the visible solar surface. 

The photosphere lies at the top of the convec- 
tion zone. This is the part of the Sun which can 
be observed in visible wavelengths, since the op- 
tical depth has approached unity in the visible 
continuum. The majority of the Sun’s energy 
that is radiated into space comes from this re- 
gion. The motions of the convection zone can 
be detected in the photosphere in various oscilla- 
tions and in turbulent granulation. Magnetic flux 
tubes appear to emerge in this region. These flux 
tubes do not emerge uniformly and are associated 
with different features seen on the Sun, such as 
active regions, coronal holes, and regions of the 
quiet Sun. The temperature continues to drop 
to a mean value of about 4500 K at the upper- 
most part of this region with an atomic density 
of 10 16 cm” 3 or a mass density of 10~ 8 gm cm' 3 . 
However, observations of IR molecular lines of 
CO suggest that there may be a significant com- 
ponent of the upper photosphere at temperatures 
no greater than 3500 K [5]. 

The chromosphere starts at the temperature 
minimum. The temperature in the chromosphere 
increases with radius to a few times 10 4 A, caus- 
ing the region to produce the hydrogen Lyman 
alpha line emission in the uppermost part. This 
temperature rise requires that mechanical energy 
be transported to and deposited in the chromo- 
sphere, since the chromosphere cannot be heated 
from below by either radiation or thermal con- 
duction. By definition, the chromosphere is hot- 
ter than most of the photosphere. The ultimate 
energy source for non-radiative heating in this 
region and the regions above is the convective 
motions in the photosphere [27, 12, 44]. The 
method by which this energy is transmitted and 
dissipated differs from region to region and is still 
being studied. The particle density varies from 
10 16 to 10 lo cm -3 in the chromosphere. 

The outermost region of the Sun’s atmosphere 
is the corona, where the typical temperature 
is observed to vary from one to several times 
10 6 K , and the electron density is between 10 7 
to 10 lo cm -3 in nonflaring regions. The coronal 
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plasma is made up almost entirely of electrons 
and ions, as essentially all neutral atoms have 
been ionized. This region is probably heated by 
either waves propagating up the magnetic field 
lines, by local current dissipation, or by some 
combination of these two mechanisms, where 
plasma motions in the photosphere are the ul- 
timate power source [57]. The magnetic field 
dominates the solar plasma at these heights with 
the plasma parameter “beta” (the plasma pres- 
sure divided by magnetic pressure) less than one. 
The corona extends outward until it merges into 
the solar wind. On large scales, and perhaps on 
smaller as well, spherical symmetry breaks down 
in the corona. 

The corona and the solar wind have no sharp 
boundary. The plasma motion in the corona 
starts a mean upward flow in coronal holes and, 
to a lesser degree, over the magnetically “quiet” 
Sun. As the collective outward velocity increases 
the solar wind “begins”, and at a few solar radii 
(depending on the underlying region) it becomes 
supersonic. This wind is known to affect the 
Earth’s magnetosphere. 

Between the corona and the chromosphere is 
the transition region, where the temperature in- 
creases from a few times 10 4 K to 10 6 K in the 
short distance of about 1000 km or less according 
to most models. There is evidence that the tran- 
sition region temperature gradient is much less in 
coronal holes than elsewhere. The transition re- 
gion is thought to be heated by energy conducted 
downwards from the more thermally conductive 
corona through the transition region, in which 
the thermal conductivity decreases downward be- 
cause of the decrease in temperature. The lower 
temperatures cause the electron thermal conduc- 
tivity to drop, since the mean-free-path of the 
cooled electrons is less. The effect of density on 
the thermal conductivity is small since density 
affects both the number of electrons conducting 
energy and the mean-free-path (See Appendix A 
for solar plasma parameters). 


1.2 Features 

The Sun’s atmosphere is not homogeneous, but 
shows distinct features which manifest them- 
selves in different ways in different atmospheric 
regions (Figure 1.2). Some of these are coro- 
nal holes, active regions, quiet regions, and polar 
plumes. The properties of many of these features 
are thought to be determined by magnetic fields. 
For further details than appear here, see again 
Stix [52]. 

Coronal holes are regions of unipolar, compar- 
atively open magnetic flux tubes, which extend 
outward away from the Sun. They do not appear 
different from surrounding areas at photospheric 
and chromospheric heights. In the corona, how- 
ever, they are dark in soft X-ray emission, be- 
cause of lower temperatures and densities than 
other solar regions. Coronal holes exist at both 
poles and often extend north or south in bands. 
Coronal holes are known to be the source of high 
speed solar wind streams. 

Active regions are regions with strong mag- 
netic fields and closed loop flux tubes, which 
reconnect to another part of the solar surface. 
They are associated with photospheric sunspots, 
which are slightly cooler than the rest of the 
photosphere and emit less visible light. They 
extend into the the chromosphere and corona 
where they are brighter than the surrounding ar- 
eas (hence have larger radiation losses). Tran- 
sient phenomenon such as flares, which release 
great amounts of energy stored up in magnetic 
fields, occur in active regions. These regions con- 
tain the hottest coronal plasma. 

The part of the Sun not covered by active re- 
gions or coronal holes is called the quiet Sun. 
Quiet Sun regions cover the majority of the Sun, 
and lie between the two extremes of active re- 
gions and coronal holes. However, there is an 
almost continuous range of smaller scale “activ- 
ity” in the chromospheric network as one goes 
from very quiet Sun to the boundaries of active 
regions. 
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Figure 1.2. A Yohkoh X-ray image of the Sun taken on October 25, 1991. The dark regions at the 
poles and extending towards the equator are coronal holes. The bright region in the lower left-hand 
quadrant of the image is a bright active region while lesser active regions are scattered across the 
solar disk. The regions of intermediate brightness are quiet Sun. The quiet Sun shows structure. 
Loops are apparent off the solar limb. 

The solar X-ray images are from the Yohkoh mission of ISAS, Japan. The Soft X-ray Telescope 
(SXT) was prepared by the Lockheed Palo Alto Research Laboratory, the National Astronomical 
Observatory of Japan, and the University of Tokyo with support of NASA and ISAS. 
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1.3 Solar Emission 

The study of emission from the Sun includes ob- 
servations over a wide range of the electromag- 
netic spectrum from radio waves to gamma rays, 
as well as measurements of solar wind particles. 
Because of the variations of emissivity and opac- 
ity with wavelength and temperature, different 
wavelength ranges allow observation of different 
parts of the solar atmosphere. The higher energy 
part of the electromagnetic spectrum can only be 
observed from space. Solar emission and many 
of the associated atomic processes are treated in 
detail in Zirin [63]. 

The Sun’s photosphere emits most of the Sun’s 
energy, and most of this emission is in the visi- 
ble part of the electromagnetic spectrum. The 
visible spectrum is made up of a continuum and 
many absorption lines. The photosphere can also 
be observed and studied in the infrared, which 
is the way we observe the deepest depths of the 
photosphere. Magnetograms are made in the vis- 
ible and infrared part of the spectrum to deter- 
mine the magnetic fields of the solar photosphere. 

The solar chromosphere was discovered by ob- 
serving the Sun during solar eclipses, which oc- 
cur when the solar disk is blocked by the Moon. 
The corona above the solar limb can be observed 
during solar eclipses as well. Coronagraphs are 
also used, which simulate a solar eclipse by block- 
ing light from the solar disk using an occulting 
disk. The Sun’s ultraviolet light output comes 
mainly from the chromosphere, with some ultra- 
violet emission coming from the transition region 
and the corona. Since this output is emission 
line dominated, the relative strengths of appro- 
priately chosen lines can be used for density and 
temperature diagnostics, as will be discussed in 
Section 2.5. 

When observed in the extreme-ultraviolet 
(EUV) the Sun emits a large number of optically 
thin spectral lines, which are produced at tem- 
peratures between several 10 4 I< to 10 6 K. Model 
calculations show that these lines are emitted 
in the upper chromosphere, the transition re- 


gion, and the corona, although most are formed 
in the corona [52]. According to these models, 
each spectral line is preferentially produced at a 
different range of temperatures, so observations 
of many lines reveal information about different 
temperature regions of the atmosphere. Analy- 
sis of these spectral lines can be used to obtain 
the density and the temperature of the emitting 
plasma. 

The corona emit in the radio range of the elec- 
tromagnetic spectrum, and the corona is opti- 
cally thick to radio waves. The background com- 
ponent of the quiet sun radiation comes from 
thermal emission [32], As the wavelength in- 
creases, the height observed increases, since the 
plasma cutoff frequency increases with an in- 
crease in density. Thus, the cutoff wavelength 
decreases with an increase in density, so shorter 
wavelengths permit us to observe closer to the 
Sun where the density is higher. When EUV 
observations are compared to radio observations, 
additional insight is gained, since both are pro- 
duced by the same plasma. For more discussion 
see Kundu [32]. 

The X-rays observed from the Sun come from 
the hotter regions of the solar corona around 
active regions and from flares. They can pro- 
vide useful diagnostics of the hot solar corona to 
determine the temperature and structure of the 
emitting region. These measurements can also be 
compared to observations of other wavelengths to 
obtain a better understanding of active regions. 

Gamma rays and hard X-rays are produced 
during solar flares and can be used as a diagnostic 
to study flares. These studies require either high 
time resolution or relatively high time and spatial 
resolution, in order to study solar flare develop- 
ment. Results obtained from these observations 
then can be used to constrain flare models. 

Measurements of the solar wind from space- 
craft which measure the mass, charge, and en- 
ergy of detected particles provide much useful 
information. These observed elemental, isotopic, 
and ionization stage abundances of the solar wind 
provide constraints on models of solar wind ac- 
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celeration, because they can be used to indicate 
coronal abundances of the source regions for the 
wind. 

1.4 Solar Corona 

Two of the fields of study of the solar corona are 
heating and relative elemental abundance, which 
are the subjects of this dissertation. These two 
subjects will be introduced in this section. The 
general coronal magnetic field structure will also 
be discussed. 

1.4.1 Heating 

Since the corona is hotter than the chromosphere 
and transition region, heat is conducted down- 
ward towards the Sun’s surface. Since energy is 
also radiated away by the solar corona, mechan- 
ical or magnetic field energy must be dissipated 
there in order to balance the radiative and con- 
ductive losses. The ultimate source of this energy 
is probably the convective motions in the pho- 
tosphere, but the method by which this energy 
is transmitted into the corona and then dissi- 
pated has not been determined. The two primary 
models are (1) the dissipation of waves propa- 
gated along the coronal magnetic field lines, or 
(2) the dissipation of currents conducted along 
these field lines. A variety of hybrid models also 
exist. See Zirker [62] for a recent survey of coro- 
nal heating theories, and the conditions a suc- 
cessful theory must satisfy. 

The coronal structure is controlled by the mag- 
netic field topology with the plasma density and 
temperature determined by the amount of energy 
deposited in the structure. If, in a given coro- 
nal structure, additional energy is deposited, the 
temperature will rise causing the temperature 
gradient to increase, which will lead in turn to an 
increase in the conducted heat flux to the tran- 
sition region. This will cause additional plasma 
to be heated and rise into the corona, increasing 
the coronal density in the structure. When less 


energy is deposited into a coronal structure, the 
plasma will cool radiatively and will sink, thus re- 
moving plasma from the coronal structure whose 
density will then decrease. 

Where the energy is deposited determines the 
position of the temperature maximum in the so- 
lar corona, and, as discussed earlier, the coronal 
density and temperature gradient. Constraints 
on coronal heating from observed densities and 
temperatures are determined from SERTS data 
in Section 5.2. 

1.4.2 Elemental Abundance Vari- 
ation 

The following brief description of relative elemen- 
tal abundance is discussed in more detail in Sec- 
tion 2.5.4 and by Meyer [39]. During the last few 
years the elemental composition of the corona has 
been determined to differ from that of the photo- 
sphere. This variation seems to be in the direc- 
tion of an enhancement of elements with low first 
ionization potential (FIP) compared to elements 
with high FIP. High FIP elements are found to- 
ward the right and top of the periodic table. This 
abundance variation is probably due to the fact 
that at temperatures around 10 4 K the low FIP 
elements are ionized, while most high FIP ele- 
ments are neutral atoms. 

Recent evidence indicates a variation in the en- 
hancement factor depending on the solar feature, 
with the highest enhancement factor of about 
10 in polar plumes found in polar coronal holes 
[59, 60]. This would indicate that the mecha- 
nism (or mechanisms) causing the enhancement 
depend on local conditions of the coronal plasma 
source sites to bring material into the corona. 
These variations could then be used as an indi- 
cator of underlying processes in the lower chro- 
mosphere, where the temperature is low enough 
for the elements with high FIP to be neutral. 

These variations show that care must be taken 
with coronal diagnostics that assume fixed ele- 
mental abundance ratios, since elemental abun- 
dance can vary. Coronal diagnostics using the 
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same element are more reliable for this reason. 

1.4.3 Magnetic Field Structures 

The magnetic field dominates the coronal struc- 
ture for the following reasons. First, as men- 
tioned earlier, energy is transmitted by waves or 
currents along the magnetic field and deposited 
into the corona. Second, electron heat conduc- 
tion dominates the conductive flux in the corona, 
therefore heat is preferentially conducted along 
the magnetic field lines. Third, the plasma is 
frozen in and travels along the magnetic field 
lines, so that each set of flux tubes is fed plasma 
independently, usually from lower regions of the 
Sun. Fourth, the magnetic pressure is usually 
much larger than the plasma pressure in the 
corona, so that different plasma densities in adja- 
cent flux tubes are possible. For a more in-depth 
discussion than is given here, see Giovanelli [24]. 

The magnetic flux tubes come in two types: 
open and closed flux tubes. Some of these flux 
tubes return to the surface of the Sun in closed 
loops and are called closed flux tubes. So called 
open flux tubes exist when the magnetic field ex- 
tends outward into the solar wind. These flux 
tubes are dragged outward by the wind. 

Some flux tubes may have a coronal plasma 
density much lower than average, or at the wrong 
temperature to emit in the observed wavelength 
band. This situation could occur if energy is not 
being deposited into a flux tube so the coronal 
plasma cools radiatively and drains out of the 
tube, or if the coronal plasma temperature in- 
creases too much. 

While the magnetic field dominates the corona 
on a small size scale, most of these variations ex- 
ist on scales smaller than the scale size of the 
SERTS observations. This causes the corona to 
appear to vary smoothly in our results, with only 
one large-scale loop observed, which will be dis- 
cussed in later sections. This has given rise to the 
concept of filling factor, which is the fraction of 
the volume filled with observed coronal plasma. 


1.5 Outline of Dissertation 

This dissertation consists of six chapters. Chap- 
ter 1 is a general introduction to the Sun, and to 
the solar corona in particular. Next, the the- 
oretical tools from solar physics are reviewed, 
which will be used to convert the observed rel- 
ative spectral line intensities into derived coro- 
nal properties. The third chapter describes the 
SERTS instrument and calibration of the instru- 
ment. Next comes the analysis chapter where 
the observed coronal density, temperature, emis- 
sion measure and elemental abundance are deter- 
mined, using the tools developed in the second 
chapter. In Chapter 5 the filling factor and heat- 
ing constraints are determined from the coronal 
properties determined in Chapter 4. The as- 
sumption of a Maxwellian electron distribution 
is examined in Chapter 6. Finally, the conclu- 
sions are given in the Chapter 7. 
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Chapter 2 
Theory 


This chapter will review the theoretical tools 
used to analyze the data and develop the coronal 
models used and discussed in later chapters. In 
the first section we will discuss the density scale 
height and how this can be related to tempera- 
ture under certain conditions. We also provide 
here a more in-depth discussion of the elemental 
composition than in the introductory chapter. In 
the second section the primary emission mecha- 
nism, collisional excitation of the observed lines, 
is discussed, along with the general physics in- 
volved in calculating the spectral line emission so 
that the reader can better understand how we do 
the calculation. Tabulated values are used where 
more accurate calculations have been done. In 
the third section, the equation for the observed 
spectral line intensity is determined by integrat- 
ing along the line of sight, assuming a hydrostatic 
atmosphere with a known scale height. The in- 
tensity is a function of several different param- 
eters such as the electron temperature, electron 
density, and the elemental abundance. Methods 
of obtaining these different parameters from ra- 
tios of spectral line intensities are then reviewed 
in the fourth section. This chapter then con- 
cludes with a discussion of the effect of variation 
of line emissivity on scale height. 


2.1 Density Scale Height 
and the Scale- height 
Temperature 

From SERTS observations, the emission from the 
solar corona is observed to decrease exponentially 
with height above the solar limb. This leads to 
the possibility of determining the corona den- 
sity scale height. The experimentally determined 
scale height can be related to a theoretical pa- 
rameter called “the scale-height temperature”, if 
we assume hydrostatic conditions as discussed in 
this section. 

The following derivation of a hydrostatic at- 
mosphere is done in spherical geometry. In a 
hydrostatic atmosphere the net force is zero, so 
the force due to gravity on a unit volume must 
be balanced by the change of pressure on that 
volume. Neglecting radial magnetic forces and 
assuming spherical symmetry, we can write 

pgAdr = — AdP , (2.1) 

where p is the mass density of the gas, A is the 
area, r is the radius, dP is the change in pressure, 
and g is the gravitational acceleration. For the 
Sun g is 

g = GA/ 0 /r 2 , (2.2) 

where G is the universal constant of gravity and 
M 0 is the solar mass. Assuming the ideal gas 
approximation. 


P = pkT/p, ( 2 . 3 ) 
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where p is the average atomic weight of the 
plasma, k is the Boltzmann constant, and T is 
the temperature, we get 


-ti&dr 

r 2 H s (T) 


where 


dp dT 

= + T ’ 

(2.4) 

R%kT 
“ pGM @ 

(2.5) 


is 


called the density scale height and has units of 
length, Rq is the solar radius. Using 


dT = ~r~dr, 
dr 


( 2 . 6 ) 


and the H' s radial dependence is assumed small 
over the distance of interest. For the observed 
conditions reported here (which will be discussed 
in Section 4.4), this effect would be less than 5%. 

The analysis described above assumes that the 
only change in the total pressure is due to the 
change in the gas pressure. This may be invalid 
if, for example, either the gradient of the mag- 
netic pressure or turbulence is important. Then 
the hydrostatic scale height would be larger than 
the effective scale height and we obtained 


i To r2 fr n 2 p dP 

IF, = HATofr + flJT + R-lpkT ~dr' 


(2.13) 


we can integrate from a reference radius r 0 to get The effective scale height depends weakly on r, 

dT but for the purpose of this dissertation this effect 

r ( ^0 + r ) — - ln( - ), (2.7) can be neglected. The scale height temperature 

Jr a \7f s (T) T ) r 2 p(r 0 ) ’ j s determined from the effective scale height by 


where p 0 is the density at r 0 . 

For an isothermal atmosphere ^ = 0 so, the 
expression simplifies to 


p = p{r 0 ) exp 


jz3-(L 

H S {T 0 ) Vr 0 



( 2 . 8 ) 


Equation 2.8 in the limit of r 0 going to infinity 
becomes the expression for the standard plane- 
parallel atmosphere. Guhathakurta et al. [25] 
uses the same equation. Since 


P = N e p, (2.9) 


T s = 


R 2 kH' 

fiGM & ’ 


(2.14) 


where H' s is determined from the radial depen- 
dence of the intensity. Even if the scale-height 
temperature differs from the electron tempera- 
ture, the density scale height is still useful in 
modeling the solar corona to obtain other pa- 
rameters. The validity of the various assump- 
tions to the observations will be discussed in Sec- 
tion 4.4.4. 


where N e is the electron density we have 


N e = N e (r 0 ) exp 


' z3 

H S (T 0 ) 




( 2 . 10 ) 


For the nonisothermal case, we derived the ap- 
proximate relationship 


p = Po exp 


P's 




( 2 . 11 ) 


where the effective scale height, 7/ a , is defined as 


To 


P's 


+ 


„2 dT 

_dr_ 


H a (T 0 )T " R%r 


( 2 . 12 ) 


2.2 Coronal Elemental 
Composition 

During the last decade it was discovered that the 
coronal elemental composition often differs from 
the photospheric composition [37]. This differ- 
ence has been determined from both solar wind 
particle measurements and emission line observa- 
tions. Since, in Section 4.6, the relative elemental 
abundance will be determined, some background 
on the present theories for variations in the rela- 
tive elemental abundance needs to be discussed, 
along with other sets of observations. Attempts 
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Figure 2.1. A schematic representation of the 
coronal elemental abundance divided by the pho- 
tospheric elemental abundance, illustrating the 
difference between high FIP elements and low 
FIP elements. (After Feldman 1992 ...)[17] 

to determine the nature of these differences and 
their causes continue. 

The elemental abundance variation seems to 
depend most strongly on the first ionization po- 
tential (FIP) of the element. In Figure 2.1 taken 
from Feldman [18] the coronal relative elemental 
abundance compared to the photospheric abun- 
dance is shown. A “step” function is observed 
with the low FIP elements (< lOeV) enhanced 
compared to the high FIP elements (> lleV). 
The size of the step is the enhancement fac- 
tor which is found to be in the range of 3-5 
[17, 38, 40]. The elements S and P, with FIPs in- 
termediate between the low FIP and high FIP el- 
ements, show an intermediate enhancement fac- 
tor. Very low FIP elements (< 6.1eV) can have 
an even greater enhancement factor than the 
“normal” low FIP elements [16, 46, 15]. 

The situation regarding hydrogen is still open 
to debate. Since it is the most abundant ele- 
ment on the Sun, it might act differently. The 
hydrogen FIP, and only ionization potential, is 
13.6 eV, which would make it a high FIP ele- 
ment. Since, at coronal temperatures hydrogen, 
is completely ionized, direct measurements are 


difficult. Initial work suggested hydrogen acted 
as a low FIP element [58]. More recent work has 
suggested that hydrogen probably acts as a high 
FIP element (see discussion by Meyer [39]), or at 
least somewhere in between (Fludra and Schmelz 
[19]). In this dissertation, we determined the ra- 
tio of the elements abundance to iron, since only 
low FIP elements are used and most of the spec- 
tral lines observed by SERTS come from iron. 
The plots and figures use Meyer’s 1985 [38] el- 
emental abundance, which treats hydrogen as a 
low FIP element. 

The enhancement factor of the low FIP to high 
FIP elements seems to vary with observation and 
the region of the corona. The enhancement fac- 
tor can be as large as 10 for polar plumes [59, 60] 
and as small as unity [48, 34]; in the quiet Sun 
upper atmosphere it is approximately 4. An en- 
hancement factor of approximately 4.5 is used 
by Meyer 1993 [39] to define his reference set of 
“coronal” abundances. Elemental abundances in 
active regions seem to vary between Meyer coro- 
nal and photospheric composition, according to 
Saba and Strong [48] and McKenzie and Feld- 
man [34]. Feldman also reports that in plumes 
and impulsive flares, elements with very low FIP 
show a factor of 2 increase over other low FIP 
elements [17]. 

These variations suggest that a mechanism ex- 
ists whereby the enhancement factor could vary, 
depending on the conditions described below, 
from one solar region or feature to another. An- 
other possibility is that there may be two or more 
mechanisms for obtaining coronal plasma from 
the photospheric plasma, one of them with low 
FlP/high FIP enhancement and a second which 
brings up less differentiated or undifferentiated 
plasma into the corona. 

One general mechanism to produce the FIP 
bias effect requires the temperature of the source 
of plasma to be less than 10 4 K and not be sub- 
jected to intense UV and X-ray radiation. Gen- 
erally, this condition is met in the chromosphere, 
photosphere, and possibly in spicules according 
to Meyer [39]. In this case, the low FIP elements 
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would be ionized and affected by electromagnetic 
forces while the high FIP elements would be neu- 
tral and not directly affected. Hence, partial sep- 
aration of neutral atoms from ions can occur. 
Then this ion rich plasma is preferentially trans- 
ported upwards so that the corona is enriched in 
low FIP elements. 

Photospheric abundances could be produced in 
the corona by a rapid heating event such as oc- 
curs during flares. In this case, the photospheric 
plasma is heated and raised into the corona be- 
fore the ions have time to separate from the neu- 
tral atoms, so no enrichment of the low FIP el- 
ements has taken place. This would just raise 
photospheric plasma into the corona and heat it 
to coronal temperatures. 

Photoionization could also affect the relative 
elemental abundance of some elements such as 
neon. A proposed mechanism is that X-ray ra- 
diation produced in the corona penetrates to 
the coronal plasma source site and ionizes neon 
atoms which have a larger ionization cross section 
than the oxygen atoms. McKenzie and Feldman 
[34] report that the relative elemental abundance 
of high FIP elements, oxygen and neon, could 
vary in active regions and flares. Schmelz [50] 
reported enhancement of neon relative to oxygen 
for certain flares. If photoionization is important, 
then some high FIP elements will sometimes act 
as low FIP elements. 

Sometimes, the relative elemental abundances 
of low FIP elements have also been observed to 
vary. The variation in low FIP elements Mg 
(7.646 ev) and Fe ( 7.87ev) found by McKen- 
zie and Feldman [34] in P78-1 SOLEX data was 
less than the uncertainties of measurement, but 
an upper limit of 1.5 was placed on the variation. 
On the other hand, Saba and Strong [48] found 
possible variations in Fe:Mg of about a factor 
of 3 using Solar Maximum Mission data. These 
two studies had a number of differences, and the 
meaning of these results has yet to be firmly es- 
tablished. 

In solar wind data, according to Meyer 1985 
[37], the elemental composition of solar energetic 


particles (SEP) is highly variable from event to 
event. This effect might result from some accel- 
eration mechanism of the solar wind, or it might 
be because of the mechanism responsible for the 
coronal composition. 

The results reported in this dissertation indi- 
cate that there is a FIP dependence of the en- 
hancement, but not as a simple “step” function. 
This will be shown in Section 4.6. 

2.3 The Atomic Physics of 
Coronal Emission Lines 

A variety of solar coronal properties such as 
plasma density, temperature, and relative ele- 
mental abundance can be determined if enough 
emission lines are available and the emission line 
formation processes are known. In the extreme 
ultraviolet (EUV) most of the emission lines are 
produced by collisional excitation, while others 
result from photoionization followed by recombi- 
nation or charge exchange. The latter processes 
are negligible for most of the cases considered 
here, however. 

Resonance scattering or photoexcitation oc- 
curs when a photon excites an ion which then 
rapidly decays, emitting another photon at the 
same wavelength but possibly in another direc- 
tion; this mechanism is discussed in detail by 
Gabriel [22] and Nakada [42], and can be as im- 
portant as electron collision excitation for suf- 
ficiently low densities or an especially intense 
background source. If emission is resonance scat- 
tered, the observed spectral line intensity will be 
affected. 

At the heights in the solar corona discussed 
here the corona is optically thin, so radiation 
emitted directly towards the observer will proba- 
bly not be affected by absorption and scattering. 
However, since some of our observed lines are 
more intense nearer the Sun than at the height 
of our observations, resonance scattering towards 
the observed can at least in principle be signifi- 
cant in those cases. This effect is most important 
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relative to emission from collisional excitation at 
larger radii where the electron density is low, 
since emissions because of the electron collisional 
excitation (discussed next) vary as the density 
squared, while resonance scattering is propor- 
tional to the first power of the density. 

Some Comments on Notation 

Before the derivation of the spectral line inten- 
sity from electron collisional excitation, the gen- 
eral notation used to indicate the spectral line or 
atomic energy level in question needs to be in- 
troduced. Since this dissertation compares lines 
from different elements, ionization stages, levels 
of excitation, and radial positions, the notation 
needs up to four subscripts. An example of the 
notation is given by E^e.xv, 5 , 1 , where E indicates 
the energy per photon, Fe is the element, XV 
is the ionization stage minus one, 5 denotes the 
initial level and 1 the ground state that the ion 
de-excites to. The general case of an arbitrary 
emission line produced by an arbitrary transi- 
tion in an arbitrary element would be written 
Ex,z,i,j where X is the element and Z is the ion- 
ization stage minus one. Densities are specified 
with three subscripts or less, for example with 
NFe,xv, 5 {r ) indicating the volume density of Fe 
XV in the fifth excitation level at radius r. The 
density of an element in a particular ionization 
stage, summed over all excitation levels will be 
Nx,z, where 


Nx,z = ^2 Nx,z,j ■ (2.15) 

j 

Finally, radial position is denoted in the common 
way as an argument in parenthesis. For example 
Nx(r) is the volume density of element X at ra- 
dius r. 

Electron collisional excitation occurs when an 
electron hits an ion with sufficient kinetic energy 
to move one of the ion’s electrons into an excited 
level. Later the electron will de-excite radiatively 
to either its ground level or to another lower level. 
This process will produce a photon at a specific 
energy, giving rise to an emission line. The emis- 
sion rate per unit volume, e, of the transition is 


e x,z,i,j = E x ,z,i,jNx,z,iAx,z,ij, (2.16) 

where E is the energy per photon, N is density 
of the emitter, and A is the Einstein spontaneous 
emission rate (s -1 ). 

Equation (2.16) can be rewritten to separate 
the terms which are dependent on the specific 
line observed and those which depend only on 
the physical conditions of the emitting plasma, 
as 


Ny 

€ X,Z,i,j = Qx,Z,i,}( T e , • /V e)-^-A r e , (2.17) 

where Qx,z,i,j(T e , N e ) is called the contribution 
function which is determined theoretically. It is 
strongly dependent on T e and normally weakly 
dependent on N e , and is defined as 


Qx,z,i, 3 (Te,N e ) 


iVx,z(T e ) 

N x 


Ex,z,ijAx,z,ij N x ,z,,(T e , N e ) 


N P _ 


N X ,2 


(2.18) 


where T e and N e are the electron temperature 
and density. The expression inside the second 
bracket of Equation 2.18 is the emission per ion, 
while the expression inside the first bracket is 
the ionization fraction. In the low density limit, 
the emission per ion is not dependent on N e , as 
will be shown in Equation 2.24. Often Q is just a 
function of T e for some spectral lines. The contri- 
bution functions for many lines have been calcu- 
lated by Mewe [35] and Landini and Monsignori- 
Fossi [41]. These authors usually present their 
calculations using an assumed set of elemental 
abundances, which replaces Qx,z,i,j(Nx /N e ) in 
Equation (2.17). They assume a single temper- 
ature and a Maxwellian electron distribution in 
their calculations. The Landini and Monsignori- 
Fossi [41] set of calculations give both density 
and temperature dependence. 

To understand the underlying physics in de- 
termining the contribution function, the general 
method of calculating it is given in Sections 2.3.1 
and 2.3.2. Section 2.3.1 discusses how to obtain 
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(2,19) 


the emission per ion, and Section 2.3.2 the ion- 
ization fraction. These two calculations are nor- 
mally done separately and often by different re- 
search groups, because the rates that affect the 
emission per ion are fast compared to the rates 
affecting the ionization fraction in the “coronal 
approximation”, where the two factors can be 
decoupled. 


2.3.1 Collisional Excitation of 
Coronal Ions 


The second part of the the contribution function 
of equation (2.18) is the emission per ion. This is 
dependent on the rate the excited level is filled, 
the rate at which electrons decay from the excited 
level to produce the emission line of interest, and 
the energy per photon as shown in this section. 
For more details, see Landini and Monsignori- 
Fossi [33] and Gabriel and Jordan [23]. As men- 
tioned earlier, calculations for different emission 
lines can be found in the literature [35, 33]. Since 
the excitation and de-excitation rates occur on 
time scales of less than a second, the effect of 
transportation of ions from one region to another 
is negligible. 

Any excited level can be populated by colli- 
sional excitation from lower levels, by radiative 
decays from higher levels, or by collisional de- 
excitation from higher ionization stage (charge 
exchange has negligible effects in the solar corona 
on the population of the excited level). The ex- 
cited level can be depopulated by collisional exci- 
tation, by collisional de-excitation, and by radia- 
tive decay to lower levels. If the ionization rate 
and recombination rate (which change the ion- 
ization state and will be discussed later) are slow 
compared to the radiative and collision rates, the 
rate of change of the excited level is 


dNx,z,i 

dt 




-N e N x ,z,i X] + X Nx,z,jAx,z,i,i 
i>* 


~N. x,z,i2lAx,z,i,j, 

3<i 

where the ground state is given by j = 1 , and C{j 
is the electron collision excitation rate (cm 3 s -1 ). 
The electron collision excitation rate coefficient 
is given by Kato [28] as 

8.63 x lO~ 6 exp(-Ex,z,i,j/kT e ) 

— pjr 7 

WiVTe 

( 2 . 20 ) 

where i > j, and Tj 7 effective collision 

strength, which can be calculated quantum me- 
chanically and is symmetric with respect to ex- 
change between i and j. By the law of de- 
tailed balance, an equilibrium argument, the rate 
of transitions from level i to level j because of 
the excitation must be equal to the rate of de- 
excitation of level j to level i [10, 6], 


Ujexp(-E x ,z y j/kT € )Cj,i = 

u,exp( -E x ,z,i/kT e )C it j, (2.21) 


where Wj is the statistical weight of the ground 
level, Ex,z,i and Ex,z,j are the energies of the 
two levels, and the collisional de-excitation rate 
coefficient can be written 


C 3 ,i = 


8.63 x 10" 6 _ 


( 2 . 22 ) 


where i > j. 

For a steady state plasma, where JA, ^ t | Z ' ' = 0, 
Equation (2.19) becomes 


iV e E 


j±i Nx,z, 3 
Ne 


Nx,z,i(Te,N e )= (2.23) 
Cj,i(T e ) + E J> , N x ,z,jAx,z, } ,t 
Cj(F e ) + Ej<, A x ,z,i,j 


To obtain a closed set of equations we include 
Equation (2.15). 

For a low electron density plasma, the first 
term of the numerator and the second term of the 
denominator dominate, with the primary source 
of excitation being from the ground level j = 1 
and with the ground level population being near 
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unity. Then Equation (2.24) reduces to what is 
called the coronal approximation 


N. 


X,Z,i 


NeC^Nxj 

J2j<i Ax,z,i,j 


(2.24) 


From this, the contribution function from Equa- 
tion (2.18) becomes 


Qx,z,i,j 


[ Ex,z,ijCubx,z,i,j ] 


Nx,z ~ 

N x J’ 


(2.25) 


where bx,z,i,j is the branching ratio 


>X,Z,i,j = 


l x,z,t,j 


J2k<t Ax,z, 


i k 


(2.26) 


This is the usual approximation used to calculate 
the intensity of a line. 


2.3.2 Ionization Balance 

The next step in determining the contribution 
function is evaluating the ionization fraction, 
i.e., the first expression in brackets in equations 
(2.18) or (2.25), which is the fraction of an el- 
ement in a particular ionization state. These 
calculations assume an equilibrium situation and 
are described in more detail by Arnaud and 
Rothenflug [3]. 

Since it is rare for an ion to doubly ionize or 
doubly recombine in one step under normal coro- 
nal conditions, the ion density is strongly depen- 
dent on the adjacent ion densities, with the rate 
of change being 

dNx 7 

— 7 f— = N e (Nx,z-i<lx,z-i + Ax,z+i q, x,z+i 

at 

—Nx,z<lx,z — Nx,z<xx,z), (2.27) 

where qx t z is the ionization rate coefficient which 
is the sum of the direct ionization rate coefficient 
and the autoionization rate coefficient, and ax,z 
is the sum of the radiative recombination rate co- 
efficients and dielectronic rate coefficients, both 
of which depend on electron temperature. (Ra- 
diative recombination dominates three body re- 
combination at densities less than 10 2O cm -3 and 


temperatures greater than 10 6 K [10].) As before, 
the plasma is assumed to be in steady state, so 
= 0, giving 

Nxz = N ^z-iqx,z-i ± Nx,z+i<xx,z+i' (2 28) 

qx,z - &x,z 

To close the set of equations given by expres- 
sion (2.28) for each stage of ionization of the ele- 
ment, we sum over all ion densities from the same 
element 

z+i 

Nx = Nx,z- (2.29) 

z=i 

Numerical ionization balance calculations have 
been done by Arnaud and Rothenflug [3] for 
various elements that are abundant in the solar 
corona. They discuss in more detail the calcula- 
tion of the ionization fraction. Arnaud and Ray- 
mond [2] have recently revised these calculations 
for iron. 


2.4 The Line of Sight Inten- 
sity 

The SERTS instrument, which is described in 
more detail in Section 3.1, can be used to ob- 
serve the corona off the solar limb, as shown in 
Figure 2.2. The emission lines observed are op- 
tically thin, so most of the light emitted toward 
the observer will reach the observer. Since the 
instrument looks along a line of sight and not at 
a particular volume, the total intensity observed 
is an integral of all the properly directed emis- 
sions along the line of sight. Since the emissions 
are dependent on the density squared, and since 
the density is falling off rapidly with height, the 
region closest to the Sun along a given line of 
sight will normally dominate. 

As just described, the observed spectral line in- 
tensity is the emission per unit volume integrated 
along the line of sight, 

tx.z.i.i = ^ j_ x 'X.ZAi dx 

(ergs cm _2 s _1 sr _1 ), (2.30) 
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where x is the path length along the line of sight 
and x = 0 is the point along the line of sight 
closest to the sun (Figure 2.2). In the coronal 
approximation, with the expression for intensity 
from Equation 2.17, this becomes 



1 f°° IV 

N'Q x ,z,ij(T e , N e )-j^-dx 

( ergs cm -2 s _1 sr _1 ). (2.31) 

As a first approximation to the elemental abun- 
dance variations between different coronal re- 
gions, we can assume that the elemental abun- 
dance along the line of sight is constant, then we 
have 


7 , _ Qx,Z,i l3 ( T e,N e ) N x ? 

h.z.,., ^ — jrj-„ N ' dz 

(ergs cm -2 s -1 sr -1 ). (2.32) 

The density-weighted average contribution func- 
tion is by definition 


n . Y7 1 j\r \ j— oo Q X,Z,i,j A' e dx 

tyx,Z,t,j( e e) joo jy 2 ^ x , 

we use the approximation 


(2.33) 


Qx,z„AT e N e ) « Q x ,z,.ATe, N e ). (2.34) 


Since the emission from electron collision goes as 
A^ 2 , the line-of-sight temperature is weighted by 
the Af? and is 


Figure 2.2. Illustrates the Line of Sight to the 
Sun and other geometrical quantities. 



n o mi ± 

JZoNidx ’ 


(2.35) 


and the line-of-sight density or line-ratio density 
which will be discussed in Section 2.5.1, is 


— njMMi 

e fZo A Qdx • 


(2.36) 


The line-of-sight electron density can be de- 
termined by numerical integration if we assume 
an exponential atmosphere of the form given by 
equation (2.10). If the scale height is less than 
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0.2 -Rq the line-of-sight electron density can be 
related to the density as a function of radius as 

TVe(r') « 0.7iV e (r), (2.37) 


where r' is the line-of-sight radius, or the mini- 
mum distance from the Sun’s center to any point 
along the line of sight. Coronal s^ale heights were 
less than 0.2R @ for all three observed regions. 

The equation for the intensity along the line- 
of-sight used in this dissertation is obtained 
by substituting Equation (2.34) into Equa- 
tion (2.32). Doing this yields 




Qx,z,;,j(7 e , Ne) Nx 
47t N e ^ 

(erg cm~ 2 sec~ 1 ), 


(2.38) 


where the column emission measure g is defined 


as 




NUx. 


(2.39) 


The column emission measure is a quantity that 
can be determined from intensity observations 
whenever reasonable estimates of the contribu- 
tion function can be obtained. The advantage 
of using a column emission measure is that the 
intensities of different lines are simply related by 
the ratios of their contribution functions. Thus, 
within the proper range of temperature and den- 
sity, we can calculate the expected line intensities 
of other lines. The density can be inferred from 
the column emission measure by modeling, but 
cannot be directly determined from the line-of- 
sight intensity of one line. See Section 2.5.1 for 
the method to determine the density from more 
than one line. The column emission measure, 77, 
can be rewritten, assuming spherical symmetry, 
as 


/ OO r 

, ^e( r ) ^2 _ r n dr - ( 2 ‘ 40 ) 
This then becomes 


’l(r') = Ni(r')L(r')f(r'), (2.41) 

where f(r') is the filling factor (Section 5.1) 
and L(r') is the effective line-of-sight length, 


and Ng(r') is the average of the electron den- 
sity squared at radius r', which is approximately 
twice the average line-of-sight electron density 
squared, (JV e ) 2 . The effective line-of-sight length, 
L(r'), can be determined by numerical integra- 
tion to within 1% for scale heights between 0.01 
and 0.4 Rq , as 


L(r') = (1.7064*//°' s -f 1.7134*//J' 5 )r' 1 ' 5 , (2.42) 


if we assume an exponential atmosphere with all 
units in solar radius units. 

Substituting Equation (2.41) into Equa- 
tion (2.38) we get 


Ix,z,iA r ') = 


<?X,Z,.',j(r e , Ne) Nx 
47 r N e 

N^(r')L(r')f(r'). 


(2.43) 


If N£ has the radial dependence of N e in Equa- 
tion 2.10, the intensity will be 


Ix,z,i,j(r') - 


Qx,z,i,j(T e , N e ) Nx 


N^(r 0 )L(r')f(r') exp 


47T 

-2 ( 1 1 


N e 


[//.(To) (r„ 




(2.44) 


The above derivation was done independently 
by the author. Guhathakurta et al. [25] derived 
a similar result including resonance scattering in- 
dicating Equation 2.44 can be written as 


Ix,z,i,j( r ') = c 

Ng(r 0 )L(r 0 , ()r n ' 5 exp 


Qx,z,t,j(r e , Ne) Nx 
Aw N„ 

-C/1 


H, 


(---) 

\r 0 r' / 


(2.45) 


to include the effects of resonance scattering, 
where C is a correction factor and the value of 
£ is between one and two because of the density 
squared dependence of electron collisional exci- 
tation and the density dependence of resonance 
scattering. 
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2.5 Spectroscopic Diagnos- 
tics 

With the development in the last section of the 
theoretical dependence of a spectral line inten- 
sity, we can now determine the properties of the 
emitting coronal plasma. This can be done by 
comparing different spectral lines which are sus- 
pected to be emitted by the same coronal plasma, 
and by determining what range of coronal condi- 
tions could produce the lines in the intensity ra- 
tios observed. Since there will always be a range 
of plasma conditions along any line of sight, val- 
ues determined from comparing different spectral 
line intensities would be an average of the condi- 
tions along the line of sight weighted by emission 
from each region. Most of the coronal lines we 
observed are produced by electron collisional ex- 
citation. The spectral line intensity from such a 
line is dependent on the local electron temper- 
ature and electron density and on the relative 
elemental abundance. 

Two general methods of modeling exist. One is 
the differential emission measure method, which 
usually assumes the low density limit, but allows 
for emitting plasma to exist at different temper- 
atures distributed as a slowly varying function. 
The differential emission measures can be de- 
termined by forcing the theoretical distribution 
of emission measures to reproduce the observed 
spectral line intensities. The result then provides 
an indication of how the plasma along a given line 
of sight is distributed as a function of tempera- 
ture. The uncertainty is dependent on the num- 
ber of spectral lines used and normally the dif- 
ferential emission measure is least constrained at 
the highest and lowest temperature ranges where 
few spectral lines exist. 

The other way to determine coronal plasma 
properties, which is used in this dissertation, 
needs fewer emission lines. In this method, we 
selectively compare the ratio of observed spectral 
line intensities whose theoretical ratio should de- 
pend only on one free parameter, in order to de- 


termine that parameter. Next, using that result, 
and then comparing lines with relative intensities 
depending on two parameters, we can determine 
the unknown parameters. This process is con- 
tinued until all the parameters that can be de- 
termined from the available data set have been 
determined. The general method for determining 
each parameter will be discussed in the following 
sections. The steps are the following: 

1. Determine N e Using spectral lines 

from same element 
and ionization stage 

2. Determine T e Using spectral lines 

from same element 
but different ionization 
stage and knowing N e 

3. Determine rj Using the spectral line 

(/ N^dx) intensity and knowing 

N e and T e and 
assuming some ion to 
hydrogen abundance 

4. Determine the Using spectral lines 
Relative Elemen- from different elements 
tal Abundance and knowing N e and T e 


2.5.1 Density-Sensitive Lines 

The ratio of intensities from two lines of the same 
ion are used to obtain the density measurement. 
Normally the line ratios are only density sensitive 
over two to three orders of magnitude of electron 
density and are approximately constant outside 
of this range. Each ion has its own range of den- 
sity sensitivity. As an example, in Figure 2.3 the 
density sensitivity of a ratio of two lines of Fe 
XIII is given [41]. These line ratios are normally 
temperature insensitive. 

The condition for spectral line ratios to be den- 
sity sensitive in dense plasma is for the collisions 
between electrons and excited ions to become 
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Fe XIII 359.63 to 359.84 Theoretical Line-Ratio Density 



Log(Den*ity) 

Figure 2,3. The theoretical ratio of intensities of 
Fe XIII 359.63 A and 359.84 A at Log(T e )=6.3. 
The ratio changes by an order of magnitude be- 
tween electron densities of 10 s to 10 lo cm -3 [41]. 

dominant. Then the first term in the denomi- 
nator of Equation (2.24) is nonnegligible. If two 
lines from the same ion are used, the density, 
N e , can be determined from the observed line ra- 
tios, which eliminates the emission measure inte- 
gral and relative elemental abundance in Equa- 
tion (2.38), so we have 

Ix,z,i,j _ Qx,z,ij(T e , Ne) ... 

Ix,z,k,i Qx,z,k,i{T e ,N e )' 

where k and / are the levels of the second emis- 
sion line. The temperature dependence in Equa- 
tion (2.46) is normally small for lines from similar 
energy levels, since the temperature dependence 
is the factor [exp( — (En — E kl )/ kT e )\. Reviews of 
calculations of line-ratio density sensitivity from 
different ions have been made by Keenan [29] 
and Doschek [14], N e (r ) can be determined from 
modeling of the solar corona. 

2.5.2 Electron Temperature De- 
termination 

After establishing an appropriate density to use, 
either by the above method or by assuming some 
reasonable density, the electron temperature can 
be determined. In this dissertation limits on the 


electron density will be determined in Section 4.2 
then the temperature will be determined. To ob- 
tain the electron temperature, one can take line 
ratios from 

lines either from the same ion or from different 
ionization states of the same element to elimi- 
nate concern of variations in the relative elemen- 
tal abundance. The disadvantages of using lines 
from the same ion are: 1) the weak temperature 
dependence of the collision rate, which requires 
a large energy difference between the originating 
levels to obtain a sensitive measure of temper- 
ature and 2) the fact that such diagnostic line 
pairs usually have a wide separation in wave- 
length and are then susceptible to experimental 
calibration uncertainties. 

The method using lines from two different 
stages of ionization of the same element, which 
depends on the temperature dependence of the 
ionization fraction, is easier to use experimen- 
tally, since here one can choose diagnostic lines 
that are close in wavelength and still have strong 
temperature sensitivity. For this method to be 
accurate, however, the observed plasma does 
need to be in steady state and reliable calculation 
of the ionization balance must exist. In addition, 
both lines need to be emitted from the same vol- 
ume of plasma, which will be shown to be the 
case in Section 4.2. The temperature-sensitive 
ratio is determined from Equation 2.38 

Ix,Z\ } i y j __ Q X ,Zl,ij{T ey N e ) 

Ix,Z2,kJ Qx,Z2 } kA T e, A r e) ’ 

where the density dependencies of Qx,z\,% and 
Qx } Z2,k are assumed to be small and the rela- 
tive elemental abundance and emission measure 
divides out. 

The uncertainties in the above equation come 
from two sources. First are the uncertainties in 
the measured intensities which are, for most cases 
in this dissertation, 15% for each spectral line rel- 
ative to one another. The second source is the 
theoretical uncertainty in the determination of 
Q. This theoretical uncertainty can be divided 
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into two parts: the uncertainty in the emissiv- the best fitted line-ratio temperature and used in 
ity and the uncertainty in the ionization balance, the rest of the analysis. 

These are difficult to estimate; in this disserta- 
tion the following values will be used: 30% for the 2.5.3 Emission Measure 
emissivity uncertainty, and (Z2 — Zl) O5 x50% for 

the ionization balance uncertainty and was deter- The emission measure, rj, can be determined us- 
mined from private communication with Arnaud ing Equation 2.38 which can be rewritten as 

[4], Raymond [47], and Monsignori-Fossi [41]. 4 tt I x ,z,i,j ,N X ,_i , 0 

To find the total uncertainty from instrumen- g = ^——^= ==y( ^ ) • [ZAV) 

tal and theoretical sources, where both are as- w ' J ' ,j 

sumed to be independent of each other and Gaus- The line-ratio density and the line-ratio temper- 
sian, add the squares of the different uncertain- ature need to be determined as described in Sec- 
ties and take the square root. The upper and tions 2.5.1 and 2.5.2. The contribution function 
lower limits to the line-ratio temperature are Q J s often only weakly dependent on N e , so only 
then determined by multiplying or dividing the a rough estimate is needed. Also the relative el- 
ratio of the intensities by the total uncertainty, emental abundance, is assumed in Section 

In Section 6.2, the log of the ratio of the intensi- 2.5.4 it will be determined. The best fit line- 
ties is shown to be approximately linearly related ratio temperature is used to determine the radial 
to the log of the ratio of the contribution function dependence of Q. 

over a range of temperatures. The uncertainty in The uncertainty in the emission measure can 
the line-ratio temperature can often be written in be best represented as an error factor, where the 
the form log(T ) = 6.30 ± 0.05 where the ±0.05 actual emission measure is within some error fac- 
is the uncertainty in the log(T). tor of the derived emission measure. The error 

Since the line-ratio temperature versus radius factor contains the uncertainties due to both the 
curves are noisy, a least-squares fit to the line- measured intensity and the contribution func- 
ratio temperature was used, where the radial de- tion. The uncertainty in the line intensity comes 
pendence of the temperature is assumed to take f ro m two sources: the uncertainty in the rela- 
the form: tive intensity of the line when compared to other 

. ("9 l' nes ) an ^ ^e uncertainty in the absolute cali- 

~~ ° a ^ r '' ' bration of the instrument, which is a factor of 

The fitted line-ratio temperature is used instead 2. The absolute calibration uncertainty will af- 
of the raw line-ratio temperature versus radius feet all emission measures from the same flight 
data to determine the “actual” scale height (Sec- equally and should not be used in the relative un- 
tion 4.4.4), the emission measure (Section 4.5), certainties of emission measures determined from 
and the radiative flux (Section 5.2). The uncer- different spectral lines. For the following analysis 
tainties of the fit are the measurement uncertain- the absolute calibration uncertainty will only be 
ties. These are then added to the instrumental included when the absolute value of the emission 
and the theoretical uncertainties to determine measure is important. 

the uncertainties in the line-ratio temperature. The uncertainty in Q is determined from the 
For the cases in this dissertation, the uncertainty uncertainty in the line-ratio temperature. The 
in T 0 is dominated by the theoretical uncertainty, total relative uncertainty in emission measure is 
while the uncertainty in a is dominated by the the square root of sum of the squares of the un- 
uncertainty in the fit. For each region analyzed certainties in Q and in the measured intensities, 
in this dissertation, multiple line-ratio tempera- The theoretical uncertainty in Q which comes 
tures are determined. One of them is selected as from the uncertainties in the ionization balance 
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and emissivity has not been used. It could be 
as large as a factor of 2. It would not be inde- 
pendent for the different stages of ionization and 
would probably have a temperature-dependent 
effect. It will be mentioned in the Table 4.7. 


2.5.4 Relative Elemental Abun- 
dance 

So far we have considered only spectral lines 
from one element. After the appropriate electron 
density and temperature have been determined, 
the relative elemental abundance of different ele- 
ments can be determined. This is simply done by 
asking what relative elemental abundance would 
be needed to produce the observed spectral line 
intensities with the emitting corona plasma at a 
particular temperature and density. This is done 
by using Equation (2.38) rewritten as 


Nx _ 47T Ix,Z,i,j 

N e Qx,z,i,j(T e , N e )r] 


(2.50) 


2.6 Line Emissivity Effect 
on Scale Height 

The original density scale height determined in 
Section 2.1 assumed a contribution function that 
is constant with radius. If this is not the case, 
the observed scale heights can disagree. By tak- 
ing into account the variation in the contribution 
function for different spectral lines as a function 
of height, we can determine the actual density 
scale height from the observed scale heights. Af- 
ter the effects of the contribution function are 
removed, the spectral lines should have the same 
scale height unless they are produced by different 
coronal structures. 

To obtain the scale height from the line in- 
tensity, we first determine the observed density 
scale height using Equation (2.43), where the ra- 
dial dependence is 

Ix.z.iA’-') oc Qx,z,ij(Z,K)r n S 

ezpl in ( i ~ ? )] - (2 - 52) 


By taking ratios of two lines from different ele- 
ments, we have 

Nxi _ Tri,zi,.,j Qx2,Z2,fc,f(T e , N e ) ,251') 

Nx 2 Ix 2 ,Z 2 ,k,l Qx\,Z\,i,}{Tei N e ) ’ 

with rj dropping out. 

The sources of uncertainties in the elemental 
abundance are the same as for the emission mea- 
sure (Section 2.5.3). However, some of them are 
not independent. The uncertainty in the abso- 
lute calibration of the instrument cancels out. 
The uncertainty due to the line-ratio tempera- 
ture from the two emission measures is not in- 
dependent and will be determined by taking the 
ratio of the emission measures at T> and T ±<tj. 
The theoretical uncertainty in Q discussed in 
Section 2.5.3 would become 2 1 * 5 and will be men- 
tioned in the Table 4.8. 


Here the radial dependence of the line of sight 
and electron density is used. If a constant con- 
tribution function, Q , is assumed, the observed 
scale height, H os is 


where 


-2 y / I \ 

H„ = ln {r^) +COnSt ’ 

(2.53) 

1 1 
y= r 

r a r' 

(2.54) 


and the radius to the 1.5 power comes from in- 
tegrating along the line of sight. 

If the contribution function is dependent on ra- 
dius, the observed density scale height will differ 
from the real density scale height. For exam- 
ple, if the contribution function decreases with 
radius, a lower scale height would be observed. 
This could occur with emission lines whose tem- 
perature of maximum abundance is cooler than 
the coronal temperature when the temperature 
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increases with radius. If the line ratio tempera- 
ture is known as a function of radius, this effect 
can be removed. First, finding the slope of the 
contribution function, c, with respect to y, Equa- 
tion (2.53) becomes 


-2 y 


H n 



+ cy- 


(2.55) 


For lines for which resonance scattering is im- 
portant, the derivation is more complicated; the 
result is that the factor of two in equations (2.53) 
and (2.55) becomes a number between 1 and 2, 
depending on the ratio of resonance scattering 
rate and collisional excitational rate according to 
Guhathakurta et, al. [25]. 


Chapter 3 
Instrument 


Before describing the analysis, a description of 
the SERTS instrument and the method of record- 
ing the raw data on film will be given. First, the 
instrument is described. Second, a general de- 
scription of how the data are obtained and then 
converted into absolute intensities are included. 


3.1 The SERTS Instrument 

The Solar Extreme-ultraviolet Rocket Telescope 
and Spectrograph (SERTS) instrument obtains 
both images and spectra of solar coronal regions. 
The SERTS instrument (see Figure 3.1) consists 
of a Wolter Type II telescope that focuses the 
light on an entrance aperture which has two wide 
rectangular apertures (lobes) connected by a nar- 
row slit. The front face of the entrance aperture 
is polished so that a reflected image of the re- 
gion surrounding what the instrument “sees” on 
the Sun is obtained for initial pointing determi- 
nation, using a broad band filter (3835 ± 30 A). 
A 1500 A thick filter of aluminum-silicon behind 
the entrance slit was used to prevent the visible 
and near-ultraviolet light from entering the spec- 
trograph. 

The extreme-ultraviolet (EUV) light which en- 
ters the spectrograph then reflects off of a near- 
normal incidence toroidal grating with a groove 
ruling density of 3600 grooves per mm. The spec- 
trally dispersed radiation is recorded on Eastman 
Kodak 101-07, EUV sensitive photographic film. 
This film is advanced, frame by frame, during the 
flight to produce multiple exposures. For each 


spectral line there appears an image of the region 
of the Sun viewed through the entrance aperture. 
The size and shape of this aperture are shown in 
Figure 3.2. The film records images in emission 
lines from 235 A to 450 A in first order, and from 
170 to 225 A in second order. In Figure 3.3 a grey 
scale picture of the SERTS 1991 spectrum from 
one of the pointing positions is shown. Figure 
3.4 shows the spectrum averaged over a number 
of rows of pixels in the wide-lobe region of the 
picture. 

A gold-coated grating was used in the 1989 
flight. In the 1991 flight, a different grating with 
a multilayer coating was used, which increased 
the efficiency in the range of 290 A to 350 A 
compared to the gold coated grating used on the 
first flight. The SERTS instrument is further de- 
scribed in Neupert et al. (1992) [43]. The SERTS 
1991 instrument modifications are described in 
greater detail by Thomas et al. 1991 [55]. 

The instrument was launched on a Terrier- 
boosted Black Brant rocket on May 5, 1989, a 
flight designated as SERTS 1989, and a modified 
version was flown on May 7, 1991, a flight desig- 
nated as SERTS 1991. During both flights, two 
different pointing positions were used, where the 
second pointing position was shifted along the 
long axis of the slit relative to the first point- 
ing position (see Figure 3.6 for the SERTS 1991 
pointing positions and Figure 3.5 for SERTS 
1989 pointing position). This allowed images of 
the region of the sun where spectra had previ- 
ously been taken and vice versa during the same 
flight. 
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Figure 3.1. The light enters the SERTS instrument through a Wolter Type II telescope and is 
focused on the polished entrance aperture which selects a region of the Sun to observe. The light 
then is reflected off a toroidal grating and the first order image is recorded on film [43]. 


3.2 Data Reduction 

The process of converting the exposed film to ab- 
solute solar intensity is a mult iple step process. It 
is necessary to understand this process in order 
to understand sources of uncertainties in abso- 
lute intensities and the limits they place on the 
analysis. 

The SERTS images were digitized by Thomp- 
son et al (1993) [56] with a Perkin-Elmer 1010M 
microdensitometer using a 10 by 10 fim pixel area 
on the film, individual pixels were then averaged 
into 50 by 50 fim pixel film-density units, corre- 
sponding roughly to 5 arc sec (3000 km by 3000 
km), comparable to the typical 5-10 arc sec in- 
strument resolution. During the analysis for this 
dissertation the 50 by 50 fim pixels were further 
averaged to 150 by 150 fim pixels. The film- 
density units were then converted into exposure 
units, using the film characteristic curve, which 
was determined experimentally also by Thomp- 
son et al. The uncertainty involved in these con- 
versions is approximately 10%. 

Next, the exposure units were converted to flux 
units by taking into account both instrumental 


efficiency and atmospheric extinction. The ef- 
ficiency of the SERTS optics was measured at 
the Synchrotron Ultraviolet Radiation Facility 
(SURF-II) of the National Institute of Standards 
an( f Technology. The instrumental efficiency was 
then determined from the combination of the ef- 
ficiencies of the different SERTS optics. The at- 
mospheric extinction is model dependent, and 
was based on the Mass Spectrometer and Inco- 
herent Scatter (MSIS) thermospheric model of 
Hedin (1987) [26]. 


Finally, a transformation from relative flux 
units to absolute solar irradiance was performed 
by normalizing the He II and Si XI 304 line in- 
tensity to previously observed, average quiet Sun 
values. The absolute calibration is estimated to 
have an uncertainty of no greater than a factor 
of two [54]. The following analysis will normally 
ignore this uncertainty, since it is the relative 
line strengths that are important in most cases. 
When this is not the case it will be noted. 
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Figure 3.2. The schematic on the left shows the dimension of the SERTS entrance aperture. For 
each spectral line, there appears an image of the region of the Sun viewed through the aperture. 
The narrow slit provide a standard spectrum of high spectral resolution. The upper and lower wide 
lobes are wide compared to the spectral line width and allow imaging of two different regions of 
the Sun. The schematic on the right shows the relative position of images from three spectral lines 
and the instrument wavelength range in first order. The widths of the wide lobes are 6.8A (not to 
scale in this schematic), and if two spectral lines lie closer than 6.8 A, then adjacent wide lobes will 
overlap. The analysis done in Chapter 4 uses the edges of the lower wide lobe, which for SERTS 
1989 and 1991 yielded images of regions off the solar limb, as shown in Figures 3.5 and 3.6. 
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Figure 3.3. A portion of the SERTS 1991 spectrum with the intensity of many spectral lines reduced 
to show features in the wide-lobe images [13]. The central region of the picture contains the narrow 
slit spectrum with spectral resolution of 50 mA. Additional spectral lines exist but are too faint to 
see for the brightness scale shown. The upper and lower wide lobes permit us to image regions of 
the Sun in many different spectra lines. Some overlapping of the wide lobes exists, but the brightest 
spectral lines produce easily observed edges. The data analyzed in this thesis were taken fiom the 
short wavelength edges off the solar limb. The solar limb starts at -8 arcmin on the short wavelength 
edge of the image and ends at -6 arcmin on the long wavelength edge of the image. The hot loop 
mentioned in Section 3.3.2 is apparent in the Fe XVI 335 A image ofT the solar limb but is not seen 
in the He II 304 A which is preferentially produced by cooler plasma. 
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Part of the Wide Lobe Spectrum 



Figure 3.4. A portion of the SERTS 1991 wide-lobe spectrum is shown above with the edges of 
the lobes from different spectral lines marked. The wavelength scale in the abscissa is aligned with 
the short-wavelength edge of the lobe. The spectrum from the wide-lobe region of the picture is 
different from a standard spectrum. The intensity of any pixel is the background intensity and the 
summation of intensities from any spectral lines for which a portion of the wide lobe is imaged 
on that pixel. The spectrum for any row of pixels in the wide-lobe region shows a series of steps 
due to the wide-lobe edges of the different spectral lines followed by corresponding decreases in 
the intensity 6.8 A away (equivalent to the physical width of the lobe). Any solar feature causes 
additional features to appear in the spectrum at the apparent wavelength corresponding to the 
appropriate physical location in the wide lobe. The intensity of selected spectral lines is determined 
from the increase in intensity across the edge of the image, as described in Section 4.1. In this thesis, 
only the line intensities at the lobe edges are used so variations across the lobe are unimportant. 
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ORIENTATIONS OF SPECTROGRAPH APERTURE 


A. ACTIVE REGION SPECTROHELIOGRAMS B. ACTIVE REGION SPECTRA 



Figure 3.5. The two pointing positions from the SERTS 1989 flight [43]. The observations analyzed 
are for the active region spectra (3.5b). The solar drawing is of the magnetic neutral lines (B z = 0). 
(Magnetic Neutral Line Drawing Courtesy of David Speich, NOAA) 
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Figure 3.6. The SERTS 1991 pointing position superimposed on a solar magnetogram [8]. The 
second pointing position was shifted along the main axis of the lobes to extend further off the 
solar limb. The white and black regions on the magnetogram indicate the sign of the line-of-sight 
component of the photospheric magnetic field. 
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3.3 Observations 


Next, we describe the conditions which existed 
off the solar limb during the two flights for the 
regions observed. The result comes from ground- 
based coronagraph observations [64, 65, 66], solar 
disk observations from a few days earlier [64, 65, 
66], and the SERTS lobe images from different 
spectral lines. 

3.3.1 SERTS 1989 

The SERTS 1989 flight had only one of the 
two pointing positions partially off the limb Fig- 
ure 3.5. Since only one exposure in this pointing 
position was made, this was the only exposure 
considered. Both edges of the off the Sun lobe 
were used for this study. The southernmost edge 
will be referred to as SERTS 1989 southern edge 
while the other edge is the SERTS 1989 north- 
ern edge, as shown in Figure 3.5 in some figures 
these were called the short-wavelength and long- 
wavelength edges. The solar sunspot activity is 
shown in Figure 3.7. The middle plot in Fig- 
ure 3.7 shows the sunspot map for the day of the 
flight. The plot on the left shows the sunspot 
map from a quarter of a solar rotation earlier, 
to show the general nature of the solar activity 
which is behind the limb on the day of the flight. 
The sunspot maps show that there are two active 
regions (AR 5464 and AR 5460) which lie near 
the line of sight. The plot on the right shows the 
Fe XIV Sacramento Peak coronagram; no Fe X 
coronagram was made for the flight date. This 
indicates that the two observed regions analyzed 
here should be classified as quiet Sun. 


the outer lobe was off the disk. The longest expo- 
sure in both pointing positions was used in this 
analysis and where some lines were overexposed, 
the second longest exposures were used. 

Only the southern edge of the lobe was used 
for this study for reasons which will be discussed 
later. This edge is called SERTS 1991 southern 
edge (see Figure 3.6). While no active regions 
exist at the solar limb, active regions exist in 
the foreground and others have rotated around 
the solar limb, as can be seen in the Boulder 
Sunspot Map from a quarter of a solar rotation 
earlier (see Figure 3.8). The Fe XIV coronagram 
shows larger emission than from the SERTS 1989 
flight. A hot loop in Figure 3.3 is visible near 
the southern edge of the lobe in the images from 
hotter emission lines. A number of foreground 
and background active regions exist. The coron- 
agrams taken for the day of the flight show both 
Fe X and Fe XIV emissions off the solar limb. 
These observations are classified as coming from 
a more active quiet sun region, since more emis- 
sion from the active regions is seen along the line 
of sight. 


3.3.2 SERTS 1991 

The SERTS 1991 flight had two different point- 
ing positions, each with 4 exposures of different 
duration. In the first pointing position, the lower 
(outer) lobe was completely off the solar limb 
with a part of the narrow slit off the disk. In 
the second pointing position, only a fraction of 


30 



April 28, 1989 


RAMEY StmSPOT 




Figure 3.7. The schematic on the left shows the active regions on the Sun from a quarter of a solar 
rotation before the flight [64, 65]. The center one shows the location of the active regions on the 
day of the flight. The one on the right shows a coronagram from the day of the flight. The location 
of the edges of the lobe can be seen on Figure 3.5b. 
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Chapter 4 
Analysis 


In this chapter, the analysis of the SERTS data 
is discussed using the theoretical tools developed 
in Chapter 2. First the criteria for line selection 
are explained. Next the density scale height of 
these lines is determined. Limits on the line-ratio 
densities are determined; then the line-ratio tem- 
perature is determined as a function of height. 
The line-ratio temperature is compared to the 
scale-height temperature after the effect of the 
radial dependence of the contribution function 
has been removed from the scale-height tempera- 
ture as was discussed in Section 2.6. Finally, the 
emission measure and relative elemental abun- 
dance are determined. 

4.1 Line Selection and 
Determination of 

Spectral Line Intensity 

The edges of the lobes from different spectral 
lines (Section 3.3) were used to analyze the so- 
lar corona off the solar limb. By comparing the 
edges of the lobe to an adjacent section of the 
film, the spectral line intensity as a function of 
position was determined. The process of obtain- 
ing spectral line intensity as a function of radius 
will be shown. Then the data sets from both 
flights will be analyzed. 

In order for the analysis of the SERTS off-the- 
limb data to be accurate, a line should meet the 
following conditions: 

• The line should be strong compared to the 


background in order to minimize uncertain- 
ties in the line intensity. 

• Adjacent lines must be resolved spectrally. 
When two ions emit at nearly the same 
wavelength, their emission lines are blended 
and only one spectral line is observed. 
Hence, this spectral line is a combination of 
the two emission lines which makes analysis 
of the combined radiation difficult because 
of the large uncertainties. 

• The emission mechanism must be well un- 
derstood. If the line emission mechanism is 
not well understood, or if the emission can 
be produced by multiple processes of similar 
strength (eg. electron collision, resonance 
scatter, and recombination), the analysis of 
the line will be uncertain. Such spectral 
lines were not included in this study. 

• Any overlapping images should be weak 
since additional uncertainties in the spectral 
line intensity will occur where the lobe im- 
ages overlap. 

After the spectral lines were selected, the 
line intensities were determined by the follow- 
ing method. The location of the lobe edge for 
each spectral line was determined from the dig- 
itized images (Section 3.2). Since the 50 by 50 
fim pixels were smaller than the instrument spa- 
tial resolution, they were merged into 150 by 150 
fim pixels. The 150 by 150 //m pixels were used 
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for the analysis in this dissertation. Thus, the in- 
tensity of the strip of 150 by 150 [im pixels along 
the edge of the lobe is then the summation of the 
spectral line intensity and the background. 

To determine the spectral line strength, the 
background is subtracted from each pixel . Since 
the instrument has some scatter which affects the 
pixels adjacent to the lobe as shown in Figure 4.1, 
the background strip was taken three 150 by 150 
fim pixels from the edge of the lobe. The spectral 
line intensity then is the intensity of the lobe’s 
edge strip minus that of the background strip. 

The solar coordinates for SERTS 1991 pix- 
els were determined by Brosius et.al. [8]. Fea- 
tures observed in He II 304 A by SERTS 1991 
were compared to Kitt Peak observation of He I 
10830 A. The scale and orientation of the SERTS 
images compared to the Kitt Peak fixed coor- 
dinate system were then determined by a least 
squares fitting routine. The solar coordinates for 
SERTS 1989 pixels were determined by a similar 
method. 

The relative uncertainty in the spectral line 
intensity is around 15% because of uncertainties 
in the intensity calibration, background subtrac- 
tion, and the determination of where the edge of 
the lobe starts. Except for weak lines, the vari- 
ation in the intensity of the spectral lines versus 
radius due to noise is normally less than the sys- 
tematic uncertainties. So, the 15% uncertainty 
is the systematic uncertainty in the overall inten- 
sity of the spectral line and not for an individual 
spectral pixel on the line. 

4.1.1 SERTS 1991 

The southern edge of the off the limb portion of 
the lobes on the SERTS 1991 flight was found 
to have six spectral lines which satisfied the se- 
lection criteria. These lines were A1 X 332 A, 
Cr XIII 328 A, Si XI 303 A, Fe XIII 359 A, 
Fe XIV 334 A, and Fe XVI 361 A. This data set 
is called SERTS 1991 southern edge. The north- 
ern edge had only one spectral line which met all 
requirements, the Fe XVI 335 A line. Since the 


SERTS 1991 Al X at 1 .05 Solar Radius 



4.1a: 



359 0 358 5 350 0 360 5 361 0 

Wavelength 


4.1b: 

Figure 4.1. The spectrum shows a fraction of 
the lobe for the Al X 332 A spectral line (4.1a). 
(The wavelength scale in the abscissa is aligned 
with the short- wavelength edge of the lobe.) Be- 
cause of scatter in the instrument, the back- 
ground needs to be taken from 4 to 6 pixels (50 
by 50 fim pixels) from the edge of the lobe. The 
Al X intensity is about 400 erg/cm 2 js/sr. For 
the Fe XIII 359 A line (4.1b), the lobe from a 
nearby spectral line forces the background to be 
taken from between 2 to 4 pixels from the edge 
of the lobe. From analysis of different lines, we 
estimate that the background is raised by 5% of 
the spectral line intensity. The edge of the lobe 
from the Fe XVI 361 A can be seen on the right 
at about 360.5 A. 
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northern edge had so few lines, only the southern 
edge wets analyzed. 

The Fe XIII 359 A had an additional compli- 
cation of a nearby lobe overlapping the region 
where the background normally would be mea- 
sured. Since the Fe XIII line was important 
for this analysis to obtain line-ratio temperature 
measurements, special treatment was given to 
this line. The background was determined us- 
ing 150 by 150 fim pixels which were only two 
pixels away from the lobe’s edge instead of the 
normal three pixels. Since this background is not 
the true background, because it is also affected by 
scattered light from the Fe XIII line, the intensity 
of the Fe XIII is therefore underestimated. To 
determine the correction factor a study of other 
spectral lines was done by taking the background 
closer to the lobe’s edge. As was done for Fe XIII, 
the correction factor was determined to be about 
1.1 (Figures 4.1a and 4.1b). The Fe XIII line was 
then assigned an uncertainty of 25% instead of 
the normal 15% because of the extra uncertainty 
and other weak overlapping lobes. 

In SERTS 1991, the overlapping portions of 
the lobe’s edge from the two different pointing 
positions were merged. This technique allows 
analysis of the solar corona to greater heights 
than either one of the pointing positions alone. 
Since the spectral lines intensity did not align 
smoothly because of calibration uncertainties, a 
box car smoothing of five pixels (150 by 150 /zm) 
was used. 


4.1.2 SERTS 1989 

The southern edge of the off-the-limb portion of 
the SERTS 1989 lobes was found to have five 
spectral lines which satisfied the selection crite- 
ria: Si XI 303 A, Fe XIV 274 A, Fe XIV 334 A, 
Fe XV 284 A, and Fe XVI 361 A. The north- 
ern edge had four spectral lines, which were Fe 
XIV 274 A, Fe XV 284 A, Fe XV 417 A, and Fe 
XVI 335 A. The reasons for the changes in lines 
between the two flights were: 


• The wider wavelength coverage in SERTS 
1989, 

• The increased efficiency and intensity of 
emission in SERTS 1991. 

A list of the lines used, the ion that produces 
each line, and its wavelength, oscillator strength, 
transition array, and the sources of the two dif- 
ferent sets of ionization balance calculations are 
given in Table 4.1. 

4.1.3 Sets of Ionization Balance 
Calculations 

Because the best ionization balance calculation 
for iron has not been determined, two different 
sets of calculations are used in this paper. The 
sources for the two ionization balance calcula- 
tions for the different elements are given in Ta- 
ble 4.1. The only difference between these two 
sets is which ionization balance calculation for 
iron was used. For the first set of ionization bal- 
ance calculations (IBALl), the results for iron 
from Arnaud and Raymond [2] are used. For 
the second set of ionization balance calculations 
(IBAL2), the results for iron from Arnaud and 
Rothenflug [3] are used. The contribution func- 
tions for the spectral lines in Table 4.1 are given 
in Figure 4.2. 

4.2 Scale Height 

The spectral line intensities decreased exponen- 
tially with solar radius. The density scale height 
was determined for each line in all three data 
sets assuming electron collision excitation to be 
the line-formation mechanism for all lines (Sec- 
tion 2.4). 

The scale-height temperatures determined 
from the density scale height using Equation 2.14 
are shown in Figure 4.3 and given in Table 4.2. 
The first two numbers in the last row in each ta- 
ble are the average density scale height and the 
scale-height temperature for the observed lines. 


35 


IBAL1 Contribution Functions IBAL2 Contribution Functions 




Figure 4.2. The normalized contribution functions are shown as functions of temperature for IBAL1 
(4.2a), and IBAL2 (4.2b). 


Table 4.1. Spectral Line Information 


Ion 

a(A) 

Oscillator 

Strength 0 

Transition Array* 

IBAL1 

IBAL2 

A1 X 

332.89 

0.287[61] 

2s 2 % - 2s2p ’P,° 

A+Ro 

A+Ro 

Cr XIII 

328.29 

0.902(36] 

3s 2 1 So — 3s3p 1 P° 

L+MF 

L+MF 

Si XI 

303.58 

0.264 [61] 

2s 2 1 S 0 — 2s2p l P° 

A-fRo 

A+Ro 

Fe XIII 

359.63 

0.048[2lj 

3s 2 3p 2 3 P, - 3s.3p 3 3 D° 2 

A-fRa 

A+Ro 

Fe XIV 

274.203 

0.200(21] 

3s 2 3p 2 P,° /2 - 3s3p 2 2 S 1/2 

A+Ra 

A+Ro 

Fe XIV 

334.171 

0.077(21] 

3s 2 3p 2 P°/ 2 - 3.s3p 2 2 D 3 / 2 

A+Ra 

A+Ro 

Fe XV 

284.16 

0.827(21] 

3s 2 ’S 0 — 3s3p ’Pf 

A+Ra 

A+Ro 

Fe XV 

417.258 

0.032(21] 

3s 2 ’S 0 - 3s3p 3 P° 

A+Ra 

A+Ro 

Fe XVI 

335.407 

0.270 [21] 

3s 2 S 1/2 - 3 p 2 P 3 ° /2 

A+Ra 

A+Ro 

Fe XVI 

360.798 

0.125(21] 

3s 2 S 1/2 - 3p 2 P,° /2 

A+Ra 

A+Ro 


Notes: 

a References for Oscillator Strength given in Brackets 
b Transition Arrays from Kelly [30] 

A+Ra. Arnaud and Raymond [2] ionization balance 
A+Ro. Arnaud and Rothenflug [3] ionization balance 
L+MF. Landini and Monsignori-Fossi [41] ionization balance 
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The Si XI density scale height temperature is the 
highest; this is probably caused by some contri- 
bution of resonance scattering to the total Si XI 
emission. The density scale heights from SERTS 
1989 are larger than for SERTS 1991, with the 
largest density scale heights occurring for the 
SERTS 1989 northern edge. 

The variation in density scale heights from dif- 
ferent spectral lines can be explained by a combi- 
nation of the change in the contribution function 
with radius because of a change in temperature 
and density of the emitting plasma, and, for some 
lines, resonance scattering (Section 2.6). There- 
fore, while the corona has an “actual” density 
scale height, the observed density scale height 
for each spectral line would vary depending on 
the radial dependence of the spectral line contri- 
bution function. In Section 4.4.4 the effect was 
modeled using the line-ratio temperature deter- 
mined in Section 4.4. The result determined in 
Section 4.4.4 is given in the last two columns in 
Table 4.2 for IBAL1 and IBAL2. The values in- 
dicated by asterisk in the last row of each table 
are the “actual” density scale heights which were 
used in the model. 

The density scale heights from the last two 
columns of Table 4.2 agree well with the observed 
density scale height for most spectral lines; there- 
fore, the assumption that the observed spectral 
lines are emitted in the same coronal structures 
is reasonable. 

4.3 Line-Ratio Electron 
Density 

Limits on the line-ratio density will be deter- 
mined in this section. The theory used in de- 
termining these limits was discussed in Section 
2.5.1. The lower density limit on N e can be deter- 
mined from comparing the ratio Fe XIII 359.63 A 
and Fe XIII 359.84 A intensities. The the- 
oretical intensity ratio changes by a factor of 
30 between electron densities of 10 8 cm -3 and 
10 lo cm -3 . Since the Fe XIII 359.84 A is ei- 


ther not observed or barely detected, the density 
needs to be greater than 10 9 cm -3 according to 
Figure 2.3. The absences of the detection of the 
Fe XIV 353 A line while the Fe XIV 334 A is 
observable indicate that the density needs to be 
less than 10 lo cm -3 . So for SERTS 1991 southern 
edge the electron density needs to be between 10 9 
and 10 ,o cm -3 . These limits were determined be- 
low I.IRq. The density will decrease with height 
as shown in Section 2.1, so both of these limits 
will decrease with height. 

An upper limit on the electron density for 
SERTS 1989 was determined by observing that 
the Fe XIV 264.82A and the Fe XIV 274.203A 
intensities were similar. This puts a strong up- 
per limit of 10 lo c77i -3 on the density. The density 
is probably much less than this but, due to the 
large background noise near the 264. 82A line, it 
was not possible to determine a smaller upper 
limit. 

4.4 Line- Ratio 

Temperature 

Next, the line-ratio temperature was determined 
using the method described in Section 2.5.2 for 
spectral lines from different ionization stages of 
the same element. Densities consistent with the 
limits on the electron density determined in Sec- 
tion 4.3 were used in determining the line-ratio 
temperature. The line-ratio temperature was 
found to vary with height as will be discussed 
below. 

4.4.1 SERTS 1991 

The line-ratio temperature was derived for the 
SERTS 1991 southern edge data set. The results 
are shown in Figure 4.4a for an electron density 
of 10 8 ' 9 cm -3 and IBAL1 (see Section 4.4.3). Re- 
sults using IBAL2 are shown in Figure 4.5a. A 
constant electron density was used, since the ex- 
act radial density dependence was not known. 
An electron density of 10 8 9 cm -3 was chosen as 
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Figure 4.3. Plots of the log of the intensity vs. radius to determine the scale height for SERTS 1991 
southern edge (4.3a), SERTS 1989 southern edge (4.3b) and SERTS 1989 northern edge (4.3c). The 
data from each ionization stage is marked by a different symbol. The fits used to determine the 
scale height are shown by the lines. The scale height and the scale-height temperature are given m 
a table on each figure. The systematic uncertainty in the intensity of the lines is 15% except lor 1-e 

XIII which is 25%. 
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Table 4.2. Temperature from Scale Height 
Table 4.2a: SERTS 1991 Southern Edge(1.05 to 1.19/?©) 


Ion A(A) 

Log(T) Max 

Density 

Scale-Height 

Predicated Density 



Abundance 

Scale Height 

Temperature 

Scale Height® 






IBALl d 

IBAL2 d 

Fe XVI 

361 

6.41(2.57e6)° 

0.10/?© 

1.4e6 

0.11 

0.11 

Fe XIV 

334 

6.26(1.82e6)° 

0 

cxj 

o 

t-H 

O 

1.4e6 

0.10 

0.10 

FeXIII 

359 

6.20(1.58e6) a 

0.09/?© 

1.3e6 

0.09 

0.09 

Si XI 

303 

6.20(1. 58e6) 6 

0.11/?© 

1.6e6 

0.09 

0.09 

Cr XIII 

328 

6.19(1. 55e6) c 

0.09/?© 

1.2e6 

0.09 

0.09 

A1 X 

332 

6.10(1.26e6) 6 

0.09/?© 

1.2e6 

0.09 

0.09 



[Avg 

0.10±0.01/?© 

1.3±0.1e6] 




NOTES: 

a Calculated from ionization balance from Arnaud and Raymond [2]. 
b Calculated from ionization balance from Arnaud and Rothenflug [3]. 
c Calculated from ionization balance from Landini and Monsignori- Fossi [41]. 

d. IBAL1 and IBAL2 defined in Table 4.1. 

e. A “actual” density scale height of 0.10/?© was used to determine the predicted density scale 
heights discussed in Section 4.4.4. 


Table 4.2b: SERTS 1989 Southern Edge(1.05 to 

1.15/?©) 


Ion A(A) 

Log(T) Max 
Abundance 

Density 
Scale Height 

Scale- Height 
Temperature 

Predicated Density 
Scale Height® 
IBALl'* IBAL2 d 

Fe XVI 361 

6.41(2.57e6) a 

O.16i? 0 

2.2e6 

0.16 

0.15 

Fe XV 284 

6.32(2.09e6) a 

o 

h— k 

© 

1.9e6 

0.14 

0.13 

Fe XIV 334 

6 26(1.82e6) a 

0.13/?© 

1.8e6 

0.12 

0.12 

Fe XIV 274 

6.26(1. 82e6) a 

0.13/?© 

1.8e6 

0.12 

0.12 

Si XI 303 

6.20(1.58e6) 6 

[Avg 

0.12/?© 

0.14±0.017?© 

1.7e6 

1.9±0.2e6] 

0.11 

0.11 


NOTES: 

a Calculated from ionization balance from Arnaud and Raymond [2]. 
b Calculated from ionization balance from Arnaud and Rothenflug [3]. 
c Calculated from ionization balance from Landini and Monsignori- Fossi [41]. 

d. IBAL1 and IBAL2 defined in Table 4.1. 

e. A “actual” density scale height of 0.125/?© was used to determine the predicted density scale 
heights discussed in Section 4.4.4. 
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Table 4.2c: SERTS 1989 Northern Edge(1.05 to 1.15/2©) 

Ion A(A) Log(T) Max Density Scale-Height Predicated Density 


Abundance 


Fe XVI 335 6.41(2.57e6) a 
Fe XV 284 6.32(2.09e6) a 

Fe XV 417 6.32(2.09e6) a 

Fe XIV 274 6.26(1. 82e6) a 

[Avg 


Scale Height 

Temperature 

Scale Height® 
IBALl 4 IBAL2' 

0.21/?© 

2.9e6 

0.20 

0.20 

O.21i? 0 

2.8e6 

0.16 

0.16 

O.15i? 0 

2.1e6 

0.16 

0.16 

O.16/?0 

O.18±O.O3/?0 

2.2e6 

2.5±0.4e6] 

0.13 

0.13 


NOTES: 

a Calculated from ionization balance from Arnaud and Raymond [2]. 
b Calculated from ionization balance from Arnaud and Rothenflug [3]. 
c Calculated from ionization balance from Landini and Monsignori-Fossi [41]. 

d. IBAL1 and IBAL2 defined in Table 4.1. 

e. A “actual” density scale heights of 0.16/2© for IBAL1 and 0.18/2© for IBAL2 were used to 
determine the predicted density scale heights discussed in Section 4.4.4. 


reasonable since it produces self-consistent plots 
for radial positions greater than 1.1 /2©. 

The systematic temperature uncertainties due 
to instrumental and theoretical uncertainties 
(Section 2.5.2) are shown in Figures 4.4a and 4.5a 
as vertical error bars in the upper left-hand cor- 
ner; they affect the overall normalization of the 
derived line-ratio temperature curves and do not 
denote relative uncertainties of each data point. 

The fit was done for the line-ratio tempera- 
ture for Fe XIII/Fe XVI above 1.1 /2©, because 
excess emission from the hot loop (Section 3.3.2) 
contributes to the Fe XVI emission below this 
height. The fits to Equation 2.48 for both the 
FeXIII/Fe XIV and Fe XIV/Fe XVI line-ratio 
temperatures are given in Table 4.3. The Fe 
XIV/Fe XVI temperature is considered to be 
more accurate because of additional noise due to 
background subtraction in the Fe XIII/Fe XIV 
line-ratio temperature above 1.1/2©. The fit for 
this ratio is plotted in Figures 4.4a and 4.5a. This 
fitted curve is used instead of the raw tempera- 
ture versus radius data to determine the “actual” 
scale height (Section 4.4.4), the emission measure 
(Section 4.5), and radiative flux (Section 5.2). 


4.4.2 SERTS 1989 

Line-ratio temperatures were determined for 
SERTS 1989 southern and northern edges in the 
same way as for SERTS 1991 southern edge and 
the results are plotted in Figures 4.4b, 4.4c, 4.5b, 
and 4.5c. The line-ratio temperatures increase 
with radius. 

SERTS 1989 data sets have a larger spread in 
line-ratio temperatures than SERTS 1991. This 
spread is within the uncertainties in the line-ratio 
temperature. The apparent spread in the line- 
ratio temperatures may be real (Section 4.4.3). 
The least-square fits to the line-ratio tempera- 
ture are given in Table 4.4 for the SERTS 1989 
southern edge, and in Table 4.5 for the SERTS 
1989 northern edge. 

For SERTS 1989 northern edge, the difference 
in the slope of the line-ratio temperature, a, for 
ratios involving Fe XV (Table 4.5) indicates that 
Fe XV 284 A is affected by resonance scatter- 
ing. This leads to more Fe XV 284 A emission 
than would be predicted by electron collision ex- 
citation. Thus, the line-ratio temperature deter- 
mined from the Fe XV 284 A /Fe XVI 335 A 
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Figure 4.4. Plots of the line-ratio temperature vs. radius for SERTS 1991 southern edge (4.4a), 
SERTS 1989 southern edge (4.4b) and SERTS 1989 northern edge (4.4c) using IBAL1 and line 
emission calculations by Landini and Monsignori-Fossi [41]. Systematic error bars are shown for 
each line-ratio temperature in the upper left-hand corner of the figure, with the symbol on each 
error bar indicating the respective line ratio. The major source of the uncertainty is from the 
intensity calibration and would shift the whole line-ratio temperature curve up or down. The line- 
ratio temperature for Fe XIII 359 A has two peaks near 1.1 These probably are not real, but 
are due to uncertainties in background subtraction. 
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Temperature Temperature 


,Serts 1989 Southern Edge Line-Ratio Temperature Log(Ne)= 9.0 
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Figure 4.5. Same as Figure 4.4 using IBAL2. The differences between the two plots demonstrates 
the uncertainties due to the calculation of the contribution functions. 
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ratio is decreased, and that determined from 
the Fe XIV 274 A/Fe XV 284 A ratio is in- 
creased compared to the same ratios containing 
Fe XV 417 A. The effect of resonance scattering 
on Fe XV 284 A is also discussed in Section 4.4.4. 

The line-ratio temperature was also deter- 
mined using the two Fe XIV lines from SERTS 
1989 southern edge. Since the differences in 
the line-ratio temperatures are within the line- 
ratio temperature uncertainties (Figure 4.4b), 
the results are in reasonable agreement. The 
Fe XIV 274 A spectral line would be influenced 
more by resonance scattering but line ratios in- 
volving Fe XIV 334 A have a higher temperature 
and reduced slope. Therefore, there is no evi- 
dence that resonance scattering contributed to 
the differences in line-ratio temperatures contain- 
ing Fe XIV. 

For SERTS 1989 northern edge, the fit of the 
Fe XIV 274 A to Fe XVI 335 A line-ratio tem- 
perature was chosen as the best fit line-ratio tem- 
perature because of smaller uncertainties in T 0 . 
For the southern edge, the average of the Fe XIV 
334 A to Fe XVI 361 A and the Fe XIV 274 A 
to Fe XVI 361 A line-ratio temperature was used 
as the best fit line-ratio temperature. The uncer- 
tainties in T 0 and a were similar except for the 
case of a determined using IBAL2 (Table 4.4). In 
that case an uncertainty of 3 x 10 5 K was used. 

4.4.3 Comparison of Different 
Ionization Balance 

Calculations 

Line-ratio temperatures determined by using ei- 
ther IBAL1 or IBAL2 are similar but do have sys- 
tematic differences. For SERTS 1991 southern 
edge (Table 4.3), the line-ratio temperatures de- 
termined by using IBAL2 (compared to IBAL1) 
were on the average about 30,000K cooler. Also, 
it was found that the slope of the line-ratio tem- 
perature versus radius was smaller by 10-20%. 

For SERTS 1989 data sets, the results are sim- 
ilar to that of the SERTS 1991. Except for 


the Fe XV/Fe XVI line-ratio temperature which 
was not available in the SERTS 1991 data set. 
The Fe XV/Fe XVI line-ratio temperature de- 
creased more than other line-ratio temperatures 
when IBAL2 was used instead of IBAL1 (Ta- 
bles 4.4 and 4.5). The result is a reduction in 
the spread of line-ratio temperatures when us- 
ing IBAL2. Therefore, the line-ratio tempera- 
tures determined using IBAL2 are more self con- 
sistent, This would indicate that Arnaud and 
Rothenflug [3] ionization balance calculation for 
iron (IBAL2) is more accurate than the newer 
Arnaud and Raymond [2] ionization balance cal- 
culation (IBAL1). However, the uncertainties 
due to the possibility of temperature variation 
along the line of sight and also in the measure- 
ments are large enough so that IB A LI can not 
be ruled out. 

4.4.4 Comparison of Line-Ratio 
Temperatures and Scale 
Height 

The variation in the density scale heights for dif- 
ferent spectral lines (Table 4.2) can be reduced 
by removing the effect of the radial dependence of 
the line-ratio temperature. The next to the last 
column in Table 4.2 contains the predicted den- 
sity scale height for each of the observed spectral 
lines using IBAL1, while the last column comes 
from IBAL2. The last row in each column is the 
“actual” density scale height which was used in 
modeling the predicted density scale heights. 

The modeling was done as follows. The model 
first uses the least-squares fit of the line-ratio 
temperatures from Tables 4.3, 4.4, and 4.5 (pa- 
rameters T 0 and a, defined in Equation 2.48) to 
calculate the contribution function, Q, of the 
different spectral lines as a function of radius. 
The observed density scale height from the spec- 
tral line for which the contribution function was 
nearly constant with radius was used as an ini- 
tial guess to the “actual” density scale height. 
This “guess” was then used to predict the appar- 
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ent density scale height for the different spectral 
lines. The initial guess was adjusted until the 
differences between the apparent and observed 
density scale heights were minimized after al- 
lowances were made for the effect of resonance 
scattering (which acts to increase the apparent 
density scale height) on some lines. 

The SERTS 1989 southern edge gave the best 
agreement, where differences between the ob- 
served and apparent density scale height were no 
greater than 0.01 Rq. The SERTS 1991 south- 
ern edge yielded similar agreement, except for Si 
XI, where the observed density scale height was 
0.02 R & larger than was predicted. The SERTS 
1989 northern edge had a larger observed density 
scale height than predicted from both Fe XV 284 
and Fe XIV 274. The Fe XV 284, the Si XI 303, 
and to a lesser extent the Fe XIV 274 A lines 
are the most susceptible to resonance scattering. 
“Actual” density scale heights from calculations 
using both IBAL1 and IBAL2 were the same ex- 
cept for the SERTS 1989 northern edge. The 
latter had the highest line-ratio temperatures, for 
which the differences between IBAL1 and IBAL2 
were greatest. 

For each observed region using IBAL1 and 
IBAL2, Table 4.6 gives the “actual” density scale 
height, the scale-height temperature, and the 
line- ratio temperature at 1.1 Rq. In these re- 
sults the scale-height temperatures and the line- 
ratio temperatures are similar, but they can dif- 
fer from each other by at least half a million 
degrees. In two of the cases the line-ratio tem- 
perature is higher, while in the third the scale- 
height temperature is larger. The primary rea- 
son for these differences is probably due to the 
assumption of spherical symmetry which is not 
completely correct. Secondary reasons for these 
differences can be turbulence, non-hydrostatic 
conditions such as flows, and divergence of the 
magnetic field, all of which are nonthermal ef- 
fects which can change the “actual” density scale 
height. It should be pointed out that, except for 
determining the spectral line contribution func- 
tions, that the two temperatures were indepen- 


dently determined. 

The results from this section indicate that the 
best lines for obtaining density scale heights are 
spectral lines that have a contribution function 
with at most a weak temperature dependence in 
the range of interest, and which are not strongly 
affected by resonance scattering. Then radial 
variation of the coronal temperature will have 
minimal effect on the inferred scale height. The 
scale-height temperatures from such lines, for ex- 
ample Fe XIII and Fe XIV in these data sets, can 
be used to give a reasonable estimate of the coro- 
nal temperature. 

4.5 Emission Measure 

The emission measure defined earlier in Equa- 
tion (2.39) can now be determined for each ob- 
served spectral line in the three data sets as de- 
scribed in Section 2.5.3. The contribution func- 
tion was evaluated using the least-squared fit 
line-ratio temperature given in Table 4.3, using 
IB ALL The error factor of the emission measure 
is given for each spectral line and is not depen- 
dent on which elemental abundance is assumed. 
The effects of using IBAL2 will also be discussed. 

The emission measures are plotted in Fig- 
ure 4.6. Meyer’s 1985 set of elemental abun- 
dances [38] was used except for the chromium 
elemental abundances where Stone’s elemental 
abundance [53] was used, since Meyer does not 
provide the chromium abundance. The emission 
measures at 1.05 Rq are given in Table 4.7 for 
four different sets of elemental abundances. The 
main difference is that the emission measures 
determined from Meyer’s 1992 set of elemental 
abundances [40] are 3.5 times smaller than for 
the other three sets of elemental abundances. 

For the SERTS 1991 southern edge the emis- 
sion measures determined for Fe XIII and Fe XIV 
(Figure 4.6a) are in good agreement. The emis- 
sion measure from Fe XVI above 1.08 Rq was 
also in good agreement with the Fe XIII and 
Fe XIV emission measures, but it was larger bc- 
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Solar Radius 

Figure 4.6. The emission measures plotted as a function of radius are shown for SERTS 1991 
southern edge (4.6a), SERTS 1989 southern edge (4.6b), and SERTS 1991 northern edge (4.6c). 
This plot uses the line-ratio temperature determined from Figure 4.4 using IB ALL The systematic 
error bars for each emission measure are shown in the upper left hand corner with a symbol in the 
center corresponding to the appropriate emission line. In Figure 4.6, the emission measure for Fe 
XVI is larger than for the other iron lines due to the hot loop. The dips in the Fe XIII emission 
measure near 1.09 and 1.12i?@ are probably due to the uncertainties in the background subtraction. 
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low 1.08 /2@. This increase in emission measure 
of Fe XVI close to the Sun is due to the hot loop 
(Section 3.3.2) which emits weakly in the other 
“cooler” spectral lines. The emission measures 
from spectral lines from other elements show the 
same general radial dependence as iron, but have 
different overall normalizations. This result will 
be discussed in greater detail in Section 4.6. 

Emission measures for SERTS 1989 are plot- 
ted in Figures 4.6b and 4.6c. The spread of 
the iron emission measures is less than the rela- 
tive calibration uncertainties shown in the figure. 
The silicon emission measure from SERTS 1989 
southern edge agrees well with the iron emission 
measures. 

Calculations using IBAL2 produced an iron 
emission measure 30% larger than for IB ALL 
Since the line-ratio temperatures determined 
from IBAL1 are systematically higher than for 
IBAL2, the emission measures of the other el- 
ements decreased by 10% for a total change of 
roughly 40% for IBALl compared to IBAL2. 

4.6 Relative Elemental 
Abundances 

The most probable explanation for the SERTS 
1991 southern edge having different emission 
measures for different elements was that the ob- 
served coronal elemental abundances were differ- 
ent from those given by Meyer 1985 [38] (for Si 
and Al) and Stone [53] (for Cr), which were used 
in deriving the emission measures. If we divide 
all the emission measures as a function of radius 
by the Fe XIV 334 emission measure we have Fig- 
ure 4.8a (IBALl) and Figure 4.9a (IBAL2) (see 
Section 2.5.4 for method and discussion of un- 
certainties). In this representation, differences 
from Meyer 1985 set of elemental abundances 
show up as offsets from a horizontal straight line 
at y=1.0. For SERTS 1991 southern edge (Fig- 
ure 4.8a and 4.9a) the Cr XIII and the Al X 
elemental abundances are high, while the sili- 
con relative abundance is low. The chromium 


and aluminum relative elemental abundances are 
roughly constant with radius, while the silicon el- 
emental abundance appears to increase with ra- 
dius. This later apparent increase is probably 
due to the increasing importance of resonance 
scattering with height for the Si XI line (Sec- 
tion 4.4.4), rather than a real change of relative 
elemental abundance with height. The observed 
relative elemental abundances are given in Ta- 
ble 4.8 for both IBALl or IBAL2, and compared 
with other published abundances for Si, Fe, Cr, 
and Al. 

One possible explanation for the apparent 
overabundance of aluminum in SERTS 1991 was 
a large amount of cooler (^ 1 x 10 6 K) plasma 
along the line of sight. One problem with this 
explanation is that cool plasma would also cause 
excess emission of Cr XIII, Si XI, and Fe XIII 
spectral lines, which is not seen. Also, the scale 
height of the cool, unobserved plasma, would 
need to be the same as for the observed plasma. 
Some rough estimates to explain just the alu- 
minum abundance would require an emission 
measure of approximately the same amount at 
1 x 10 6 K as at 2 x 10 6 7T, which can be ruled 
out based on the low flux in the other low- 
temperature spectral lines. 

For the SERTS 1989 flight, only the relative 
abundance of silicon-to-iron can be determined 
because the lines from other elements were too 
weak to observe. The results are given in Fig- 
ure 4.8b for IBALl and Figure 4.9b for IBAL2. 
The silicon-to-iron ratio does agree with Meyer 
1985 [38], using the same emission lines for which 
SERTS 1991 implied a low silicon-to-iron ratio. 

An estimate of the aluminum to iron relative 
abundance can be made for SERTS 1989 south- 
ern edge. In Figure 4.10, SERTS 1989 spectra at 
1.1 Rq are compared to the SERTS 1991 spec- 
tra scaled by the relative emission measures and 
taking into account the different line-ratio tem- 
peratures. The SERTS 1989 aluminum emission 
shows up at the three sigma level of the back- 
ground noise which is the solid horizontal line 
crossing the figure. The aluminum intensity is 
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Table 4.7. Emission Measure 


At l.O5/? 0 in units of 10 27 cm 5 assuming IBALl^ 


Table 4.7a: SERTS 1991 Southern Edge for T(r)=1.90e6+0.59e6*(r-l) 


Elemental 

Abundances 

Fe XVI 

Fe XIV 

Fe XIII 

Cr XIII 

Si XI 

A1 X 

Allen [1] 

10.7 

11.9 

12.7 

17.9 

7.1 

102. 

Meyer 1992 [40] 

3.2 

3.6 

3.8 

C 

1.7 

21.1 

Meyer 1985 [38] 

10.8 

12.0 

12.9 

c 

6.0 

91.5 

Stone [53] 

10.7 

11.9 

12.7 

22.1 

7.5 

95.7 

Error Factor 0 

1.23 

1.15 

1.26 

1.17 

1.18 

1.19 


NOTES 

a. The error factor (“within a factor of ...”) is a combination of the instrument uncertainty and 
the uncertainty due to the line-ratio temperature. It does not include the absolute flux 
uncertainty of a factor of 2. In addition, the ionization balance calculations and emissivity 
calculations could have a worst case uncertainty as large as a factor of 2, which is not included in 
Table 4.7. This latter uncertainty however, will not be independent of the temperature 
uncertainties. Since the uncertainties in the temperature and in the ionization balance and the 
emissivity are interconnected, we can reasonable replace the error factor in the above table by a 
factor of 2, which can then be viewed as an upper limit on the total uncertainty. 

b. The Stone set of elemental abundances [53] does not contain hydrogen, log(AV e /Ab/) was set to 
the value in Allen [1] of -4.4 for comparison. 

c. The Meyer papers [38, 40] do not have an abundance for chromium. 

d. IBALl defined in Table 4.1. 
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Table 4.7 cont’d: Emission Measure 
At 1.05/1© in units of 10 27 cm -5 assuming IBALl d 
Table 4.7b: SERTS 1989 Southern Edge for T(r)=1.84e6+0.70e6*(r-l) 


Elemental 

Abundances 

Fe XVI 361 

Fe XV 284 

Fe XIV 274 

Fe XIV 334 

Si XI 

Allen 

1.03 

0.72 

0.94 

0.72 

1.06 

Meyer 1992 

0.31 

0.22 

0.28 

0.22 

0.26 

Meyer 1985 

1.04 

0.73 

0.95 

0.73 

0.90 

Stone c 

1.03 

0.72 

0.94 

0.72 

1.12 

Error Factor* 1 

1.56 

1.31 

1.15 

1.16 

1.27 

Table 4.7c: SERTS 1989 Northern Edg< 

e for T(r)=2.02e6+1.25e6*(r 

-1) 

Elemental 

Abundances 

Fe XVI 335 Fe XV 284 Fe XV 417 

Fe XIV 274 


Allen 

0.985 

0.662 

0.586 

0.619 


Meyer 1992 

0.295 

0.198 

0.176 

0.185 


Meyer 1985 

0.997 

0.671 

0.594 

0.627 


Stone 6 

0.985 

0.662 

0.586 

0.619 


Error Factor 

° 1.63 

1.24 

1.22 

1.48 



NOTES 

a. The error factor (“within a factor of ...”) is a combination of the instrument uncertainty and 
the uncertainty due to the line-ratio temperature. It does not include the absolute flux 
uncertainty of a factor of 2. In addition, the ionization balance calculations and emissivity 
calculations could have a worst case uncertainty as large as a factor of 2, which is not included in 
Table 4.7. This latter uncertainty however, will not be independent of the temperature 
uncertainties. Since the uncertainties in the temperature and in the ionization balance and the 
emissivity are interconnected, we can reasonable replace the error factor in the above table by a 
factor of 2, which can then be viewed as an upper limit on the total uncertainty. 

b. The Stone set of elemental abundances [53] does not contain hydrogen, log( N Fe /N H ) was set to 
the value in Allen [l] of -4.4 for comparison. 

c. IBAL1 defined in Table 4.1. 
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Relative Elemental Abundance SERTS 1991 Southern Edge 



1.00 1.05 1.10 1.15 1.20 

Solar Radius 


Relative Elemental Abundance SERTS 1989 Southern Edge 



1.00 1.05 1.10 1.15 1.20 

Solar Radius 


Figure 4.8. Plots of the relative elemental abundance for SERTS 1991 southern edge (4.8a), and 
SERTS 1989 southern edge (4.8b) using IB AL1. These plots used the emission measures determined 
earlier divided by the emission measure of the Fe XIV line for Figure 4.8a and the Fe XV for Figure 
4.8b. The systematic error bars are shown in the upper left hand corner. 
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Relative Elemental Abundance Relative Elemental Abundance 


Relative Elemental Abundance SERTS 1991 Southern Edge 



Relative Elemental Abundance SERTS 1989 Southern Edge 



Figure 4.9. The same as Figure 4.8 except using IBAL2. 


4.9a: 


4.9b: 
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0.63 ± 0.21 of the scaled SERTS 1991 aluminum 
intensity. Part of the SERTS 1989 aluminum 
intensity probably comes from cooler emissions 
which can be determined by noting that the 
SERTS 1989 emission in the 320 to 325 A region 
is greater than the scaled SERTS 1991 intensity 
in the same region. The 320 to 325 A region 
has lobes from Si VIII and Mg VIII lines which 
have temperatures of maximum abundance at 
800,000K. Plasma at this temperature would also 
emit A1 X 332 A. Therefore, the SERTS 1989 
southern edge aluminum abundance is probably 
even less than the upper limit (0.63 ± 0.21 of 
the aluminum abundance determine for SERTS 
1991) determined from Figure 4.10. 

We discussed the importance of the FIP of dif- 
ferent elements and why one might expect a dif- 
ferential abundance between high and low FIP el- 
ements separated in the chromosphere, in Section 
2.2. It is interesting that, for the SERTS 1991 
southern edge, the relative elemental abundances 
increase for elements of decreasing FIP. The 
Meyer 1985 set of coronal abundances [38] has 
low-FIP elements (below < lOeV) more abun- 
dant than the high-FIP elements by a constant 
enhancement factor, when compared to photo- 
spheric abundances. All of the spectral lines an- 
alyzed come from low-FIP elements, but instead 
of a simple step function describing the abun- 
dance variations as Meyer suggested, an enhance- 
ment dependent on FIP exists for the SERTS 
1991 southern edge. However, the SERTS 1989 
southern edge set of elemental abundances indi- 
cates a constant enhancement factor. 
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100 


Estimate of the SERTS 1989 Aluminum Abundance 



Figure 4.10. The 320 A to the 335 A spectrum from SERTS 1989 (solid) and SERTS 1991 (dotted), 
where the SERTS 1991 intensity was divided by a factor of 14.0 due to the differences in the emission 
measures and the effect of the line-ratio temperatures on the aluminum contribution function. The 
A1 X 332 A and the Fe XIV 334 A are detectable in both plots with a vertical line indicating the edge 
of the lobe. The horizontal line indicates the 3-sigma level detection above the background for a 
lobe which is nine pixels (50 by 50 /xm pixels) wide. If the aluminum abundances for the two flights 
were the same, the SERTS 1989 and the SERTS 1991 intensities would be in agreement. Since 
the SERTS 1989 intensity is less, the aluminum abundance is less for SERTS 1989. The aluminum 
line can be efficiently produced by a one million degree plasma. The SERTS 1989 seems to have 
relatively more million degree plasma than SERTS 1991, since the 320 to 325 A region has relatively 
more emissions where a number of Si VIII and Mg VIII spectral-line lobes exist. Therefore, this is 
an upper limit on the SERTS 1989 aluminum abundance. 
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Table 4.8: Relative Elemental Abundances 


o ^ 50 

£ g - 


O O e 

o o 
o o 


00 

A1 O ^ 

CO O - <0 
00 5 


o> § 

5 ? § 


o M to 
01 g ^ o 
<M 0 <M 


S 0 

O o 
C- O 


10 2 ^ ^ 
^ 2 ^ ^ 
co 0 


to n h- 10 
h 00 h- O 
00 N 6 lO 


u x >> 5 ; 
^ cn .£ ^ 

* 2 I § 

c .2 $ -o 

.sals 

g * B * 

<u «3 . 

'r -C ^ 1? 

E ^ > U3 

O c; <v 

* .2 2 a 

O *3 .2 

.2 c 5 <» 
*. 0 >* 
5gS.S 

w C O nj 

»s g •■§ -e 

^ oi 5J 
& ^ u 5 

Lj -W ,»— • 3 

£ ■ < <$ 

v v 

a s* o 
H ^ -r -2 

u -ij > ns 

. 2 '« « 
tn .2 _C 

« c. 5 -*> 
| E <u £ 
-E « -n % 


c 

<y 

(m 

eti 

&* 

.2 

d 

u> 

w- 

w 

c 

.2 

CL 

s- 

*0 

O 

<L 

%-3 


jC 

c 

Jt 

0 



D 

si 

.£ 

0) 

J= 


~o 

C 



u 

it 

t- 

(t 

0 
* 1 

4) 

U 

c- 

~c 

a; 


^ 3 S 

HD ^3 ~ 
. _S y 


f— <( 

L— 


c • 

* 

O 


.2 " 

^4> 

(-< 

O 

’ 5 

4-9 

t_i 

<e 

JO 

u 

cd 

L* 

0) 

U 

c 

D 

N 

*s ■ 
.2 * 

_c 

ct3 

<U 

a; * 

j*-f *» 


o c; 
c ^ 

y) CN 
<U 

O <r- 

-o 3S 

, CQ 


C 5 H * 
5 ^ • °2 

• ~ C C-l 


E E 

3 a 

E | 
C 2 I 

O J3 — 

i o < 


tn « ‘-J 
W 

H 5 2 
o ^ 

2 «i i 


63 



Chapter 5 

Implications of the Coronal Parameters 


The coronal parameters determined in Chap- 
ter 4 have implications for determining the fill- 
ing factor, minimum electron density, and heat- 
ing requirements for the solar corona, and will be 
discussed in the following sections. 

5.1 Filling Factor and Min- 
imum Electron Density 

Since the emission measure is determined over 
a field of view times a depth, the small scale 
structures are averaged in this observed vol- 
ume. While the line-ratio density determines the 
weighted average density of the emitting struc- 
tures, the latter may not completely fill the po- 
tentially observed volume. The filling factor is 
the volume of the emitting structures divided by 
the observed volume. If the filling factor is less 
than unity, then only a fraction of the observed 
volume is filled with plasma emitting in the ob- 
served wavelength. If the filling factor is one, the 
observed volume is uniformly filled. 

The average of the electron density squared, 
/V?, is related to the emission measure by Equa- 
tion (2.41), which can be rearranged as 

(5 - 1) 

The minimum electron density necessary to ex- 
plain the observed emission measure will then be 

{N e ) min = yj±. (5.2) 


The minimum electron density at 1.05 R & is 
given in Table 5.1. In the table the emission 
measures, ?/, determined in Section 4.5, and the 
density scale heights used to calculate L were de- 
termined in Section 4.4.4, IBALl, and the Meyer 
1985 elemental abundances [38] were used. The 
minimum electron density using Meyer 1992 el- 
emental abundances [40] would be reduced by a 
factor of 0.55. 

The limits on the filling factor can be deter- 
mined by comparing the minimum electron den- 
sity and the limits on line-ratio density (Sec- 
tion 4.3). The line-ratio density is an average 
density along the line-of-sight and not the den- 
sity at a single radial position. For an exponen- 
tial atmosphere with a density scale height of 
less than 0.2 i?@, the line-ratio density times y/2 
would be the electron density at the point along 
a line-of-sight closest to the center of the Sun 
(Equation 2.37). Substituting Equation (5.2) 
into (5.1), the filling factor would then be 



where N e is the line-of-sight density determined 
from the line-ratio density (Equation 2.36). The 
limits on the filling factor are given in the last 
column of Table 5.1. 

Since the line-ratio density for SERTS 1991 be- 
low 1.1 was determined to be between 10 9 
and 10 lo cm~ 3 (Section 4.3), this implies a filling 
factor between 6 x 10” 4 to 6 x 10~ 2 , for IBALl 
and Meyer 1985 elemental abundances [38]. The 
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filling factors for both SERTS 1989 southern and 
northern edges need to be greater than 3 x 10r 5 
because of an upper limit on the line-ratio den- 
sity of 10 lo cm -3 . The filling factors for SERTS 
1989 are probably much greater than the lower 
limit, because the maximum density estimate is 
probably too high. 

5.2 Heating 

Since the corona is much hotter than the photo- 
sphere, the corona requires a source of energy to 
be deposited in it by some non-radiative mech- 
anism or mechanisms. To estimate the energy 
deposited in the corona, one has to assess the 
energy leaving the corona by heat conduction to 
the transition region, the energy radiated away 
in the electromagnetic spectrum, and that de- 
posited into the solar wind. If a steady state 
exists, the known energy leaving a region of the 
corona is equal to the energy deposited in that 
region. This estimate of the required energy de- 
position will then impose constraints on coronal 
heating models. 

The temperature gradient in the lower corona 
causes heat conduction to transport energy to 
the transition region. Electron heat conduction 
(F c ) dominates the conduction at coronal tem- 
peratures, and along the magnetic field line. The 
conductive flux, according to Spitzer [51], is 

F c & lQ~ 6 T 2 ' s -f-(ergs cm~ 2 s ~ l ) (5.4) 

dz 

for a Maxwellian particle distribution. In equa- 
tion (5.4) T is the temperature, and z is the dis- 
tance along the magnetic field line. Across the 
field lines, the heat conduction is negligible. For 
example for N e = 10 9 cm -3 , B — 10 gauss, and 
T = 2 x 10 6 K the ratio of the conductive flux 
across magnetic field lines to along the field line 
is lO” 10 . 

For the SERTS data, any change of heat con- 
duction over the height observed was less than 
the uncertainty of the measurement. Limits on 


the change in heat conduction were determined 
by the following means. The electron heat con- 
duction equation as a function of height can be 
written as a linear equation 


Fc{z) = F 0 + Z F' c . 

(5.5) 

Substituting equation (5.4) for F c 


10 - 6 r 2 .s^ fa F 0 + z FI 
dz 

(5.6) 

and integrating with respect to z 

4 


10 -6 ,F 

(5.7) 

T 3 « F a z + z 2 + const 

o.o z 


F 0 can be determined from Equation 5.7 and the 
line-ratio temperature determined in Section 4.4. 
Then an upper limit on the divergence of heat 
flux, F c \ (though zero change of heat conduction 
fits the data) can be determined. The results are 
given in Table 5.2. 

Next, the radiation per unit volume must be 
estimated from the observations. The radiative 
power, F r , can then be estimated from 

F r = e T {T)N* ( erg cm' 3 s _1 ), (5.8) 

where ej is the total radiative loss function, 
which is a function of the temperature and the 
elemental abundances [9]. Substituting equation 
(5.1) for N*, equation (5.8) becomes 

F T = e T v(r)/fL(r,H s ). (5.9) 

Unlike heat conduction, the radiation per unit 
volume of a structure depends on the filling fac- 
tor, as indicated explicitly in the last column of 
Table 5.2. 

The radiative loss function is tabulated by 
Cook [9] who used Meyer 1985 elemental abun- 
dances [38]. The radiative loss function can be 
relatively easily determined, using Cook’s paper, 
for other sets of elemental abundances. Since 
iron is the primary emitter in the 2 x 10 6 K tem- 
perature range, the value of F r will change by 
only 5-10% for the different sets of elemental 
abundances used in this dissertation. 
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Table 5.1. 

Data Set 

SERTS 1991 Southern Edge 
SERTS 1991 Southern Edge 
SERTS 1989 Northern Edge 


Minimum Electron Density and Filling Factor 
Minimum Electron Line-Ratio Filling 


Density 0 ’ 6 
3.6 ±0.6 x 10 8 

8.3 ± 1.3 x 10 7 

7.4 ± 1.2 x 10 7 


Density 

Factor"’ 6 

io 9 -io 10 

0.0006 - 0.06 

< io 10 

> 3± 1 x 10“ 5 

A 

l—* 

o 

© 

> 3± 1 x IO” 5 


NOTES: 

a: For Meyer 1985 elemental abundances [38] and IBAL1. 

b: A systematic uncertainty of a factor of 2.0 in the filling factor, and a factor of 1.4 for the minimum 
electron density exist because of the absolute calibration (Section 3.2.) 


Table 5.2. Comparison of Radiation to the Divergence of the Heat Flux 




range 

Tf 

o 

H 

X 

F' c x IO" 6 

F r x IO" 6 // 



Rq 

(erg/cm 2 j s) 

(erg /cm 3 / 's) 

(erg/cm 3 /s) 

SERTS 

1991 south side 

1. 1-1.2 

6 ± 2 

0 ± 2 

45. ±9. 

SERTS 

1989 south side 

1.02-1.15 

5 ± 1 

0±1 

1.8 ±0.4 

SERTS 

1989 north side 

1.02-1.15 

14 ±2 

0 ±2 

0.98 ± 0.28 


Note: F r determined at 1.05 R@. F a , F' c , and F r are respectively the constant part of the divergence 
of the heat conduction, divergence of the heat conduction, and the radiative power. 


Since the temperature increases with radius 
(Figure 4.4), heating above the observed height 
range is needed. The energy deposited in the so- 
lar wind is an extra energy requirement of the 
energy deposited above the observed region. 

In Figure 5.1, heating is shown to be neces- 
sary for at least part of the observed range for a 
filling factor of unity for SERTS 1989 southern 
edge and SERTS 1991 southern edge. Heating 
over the entire observed range for all three cases 
would be required for a filling factor of 0.1 or less. 
An upper limit on the filling factor of 0.06 was 
determined for SERTS 1991 southern edge (Ta- 
ble 5.1), while no upper limit for SERTS 1989 
was determined. A filling factor of 0.1 or less 
is reasonable. For more discussion of the filling 
factor see Section 5.1. 
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10 F 


Qomparison of Radiation to the Divergence of the Heat Conduction 
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Figure 5.1. Heating Constraints 

The heating constraints are shown for SERTS 1991 southern edge (5.1a), 
SERTS 1989 southern edge (5.1b), and SERTS 1989 northern edge (5.1c). 
The maximum heating of the corona by heat conduction (horizontal line) is 
compared to the radiative flux as a function of radius for a filling factor of 
unity (+) and 0.1 (asterisk). 


5.1a 
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Power(erg/cm A 3/sec) Power(erg/cm A 3/sec) 



5.1b: 


Comparison of Radiation to the Divergence of the Heat Conduction 

10 F SERTS 1'989 Northern Edge, Meyer 1985 Elemental Abundance, IB Alii 
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Chapter 6 

Upper Limits on Bi-Maxwellian 
Electrons 


A Maxwellian distribution for the electrons ve- 
locities has been assumed in the previous chap- 
ters of this dissertation since, given enough time, 
the electrons will reach a Maxwellian distribu- 
tion in an isolated plasma. In this section we 
will determine if this is a good assumption for 
our analysis. 

Bi-Maxwellian distributions can be produced if 
high-energy electrons, produced in one region of 
the atmosphere, travel to another region without 
undergoing collisions. This is due to their longer 
mean-free path. The mean-free path is 

Xmfp = ~Na 

where a is the cross section and N is the density 
of particles. Scattering due to many cumulative 
small-angle scatterings is dominant. According 
to Krall and Trivelpiece [31] the cross section 
for 90° scattering due to many cumulative small- 
angle scattering events is 


2.6 x 10- 12 

E 2 


-cm 


(6.2) 


where E is the particle kinetic energy in eV. So 
we get 


'm/p ~ 


3.8 x lO 11 ^ 2 

N 


-cm. 


(6.3) 


The higher the energy of the electrons, the longer 
the mean free path will be. So, the high energy 
electrons produced by flares, nanoflares, or by 


other means in one region of the corona can prop- 
agate to another region of the corona about the 
distance of one mean -free- path. In the corona 
typical nonflaring kinetic energy for electrons is 
on the order of 200e V ; these electrons would have 
a mean-free-path of 150 km for N e = 10 9 cm~ 3 . 
Electrons produced in flaring regions where the 
thermal electrons could have kinetic energies of 
1 OOOe would have a mean-free-path of 3800 
km for the same electron density this is still 
small compared to the hydrostatic scale height 
of 0.1 Rq or 64,000 km. The high energy elec- 
trons then can have a nonlocal origin compared 
to the lower energy electrons. 

This process could lead to a local bi- 
Maxwellian distribution of electrons with an ex- 
cess number of high energy electrons for the lo- 
cal “mean” temperature. An upper limit on the 
excess number of high energy electrons can be 
determined both from theory and from observa- 
tions. Limits determined from both methods are 
provided in this chapter. 

The high-temperature electrons will loose en- 
ergy by collisions on the time scale given by 
Spitzer [51] as the time of equipartition, repro- 
duced below as Equation (6.28). For coronal con- 
ditions, this cooling time is on the order of tenth 
of a second. Since an equilibrium situation is as- 
sumed and the ionization balance time scales are 
on the order of tens of minutes, the limits apply 
to a continuous injection of high-energy electrons 
and not a single, transient injection of electrons. 
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6.1 Symbols 

This chapter uses the same sets of symbols that 
are used in the rest of the dissertation except for 
the following. Qx,z is the ionization cross sec- 
tion, E is the kinetic energy of the electron, a is 
the mixing fraction, C and D are fitting param- 
eters, and Ti is the measured line-ratio tempera- 
ture. 

6.2 Observations 

The effect of an excess number of high-energy 
electrons on the spectra will be analyzed to de- 
velop a diagnostic test. When an electron collides 
with an ion, in general, one or a combination of 
four things can happen to the ion: 

1. The electron is elastically scattered with 
only kinetic energy transferred; 

2. the ion can be excited; 

3. the electron can be captured so the ioniza- 
tion stage of the ion will decrease; or 

4. the electron can ionize the ion so the ioniza- 
tion stage of the ion will increase. 

The higher the energy of the electron, the more 
likely possibility 4 will occur. So the addition of 
high-energy electrons will preferentially increase 
the ratio of a higher ionization stage compared 
to a lower ionization stage. This would have the 
effect of raising the apparent temperature of the 
plasma determined by line-ratio diagnostics (de- 
scribed in Section 2.5.2). 

To model the effect on the ratio of intensities of 
spectral lines from different ionization stages we 
will make the following assumptions. The elec- 
tron distribution will be modeled as a sum of two 
Maxwellian distributions 

FAEJ^T^a) = (1 -a)F(E,T l ) + aF(E,T 2 ), 

(6.4) 

where T x and T 2 are the temperatures of a cool 
and a hot component respectively, and a is the 



Figure 6.1. The net electron distributions for two 
Maxwellian components with kT\ = 190eV and 
kT 2 = 2000eV, respectively, for two mixing frac- 
tions (a = 0.1 is shown with a solid line and 
a = 0.01 with a dashed line, and the summations 
are shown as dotted lines). 

mixing fraction of the two components (i.e., the 
fraction of the total number of electrons in the 
hotter component) (see Figure 6.1). F(E,T l ) is 
a normal Maxwellian distribution written in en- 
ergy space as 

(65) 

A two-component electron distribution where 
each component is a Maxwellian makes the anal- 
ysis readily tractable, and so is adopted only for 
illustration and to obtain a rough estimate on the 
excess number of electrons. 

In steady state the number of ions being ion- 
ized needs to match the number of ions recom- 
bining, so neglecting double ionization or recom- 
bination (see Section 2.3.2) we have 

Nx,z<lx,z(F{E)) = Nx,z+\ a x,z+\(F{E)), 

( 6 . 6 ) 

where Nx,z is the density of ionization stage Z 
for element X , qx,z(F(E)) is the ionization rate, 
and ax,z+iF(E ) is the recombination rate. Both 
q x z and Q x,z+ 1 are dependent on the electron 
energy distribution. Equation (6.6) can be rear- 
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ranged as 


Collision Cross Sections 


Nx,z + 1 _ qxMF^E)) 
Nx,z &x,z+\{F{E)) 


The addition of an excess number of high en- | 
ergy electrons will change qx,z and atx,z+i , so the g 
ratio of densities of adjacent ionization stages will f 
change. The ionization rate (see, e.q., Arnaud § 
and Rothenflug [3]) is given by 


qx,z(F(E)) = r Qx,z{E)F(E)dE, ( 6 . 8 ) 

Jo 



where F(E) is the electron distribution and Qx,z 
is the total ionization cross section (Figure 6.2). 
Next, using the electron distribution given in 
Equation (6.4), we get 


Figure 6.2. The ionization cross section as a func- 
tion of energy shown for Fe XV, Fe XVI (solid) 
and Ni XVII(dotted). 


<Ix,z(Ti,T- 2 ,a) — 

(1 -«) I" Qx,z(E)F(E,T l )dE 

Jo 

+« [°° Qx,z(E)F(E,T 2 )dE. (6.9) 
Jo 

Each integral can now be replaced using Equa- 
tion (6.8) giving 


A similar relationship for the recombination 
rate can be derived , giving 


ftx,z+i(Ti,T2,«) 

ox,z+i(Fi) 


1 + u~fa(Ti, T 2 ), 


(6.14) 


qx,z{T\, 7), a ) = (1 — a)qx.z{T 1 ) -f a qx,zO 2 ), 

( 610 ) 

where qx,z(T\) and qx.zi'l 2) are the ionization 
rates over a single Maxwellian distribution and 
can be determined, for example, by using Arnaud 
and Rothenflug [3] or Arnaud and Raymond [2] 
calculations and qx,z{T\, T<i, a) is the ionization 
rate for the electron distribution given in Equa- 
tion (6.4). 

Equation (6.10) can be rearranged to give 


gx,z(Fi,T 2l Q) , / qxAJA _ \ 

«x.z(Ti) WATi) ) 


( 6 . 11 ) 


Next, 7^ is defined as 

qx,z{U) 

So, Equation (6.11) becomes 

<1X,z{F\iT2iO.) 


qxA'i 1) 


( 6 . 12 ) 


= 1 + aigiTu T 2 ). (6.13) 


where 7 a is defined as 

7»(7’.,T!) = 2i 4^ 7-1. (6.15) 

a X,z(j l) 

In both cases the mixing ratio of the two compo- 
nents has been separated from the temperature 
dependence. 

The ionization rate is primarily determined by 
the fraction of the electrons with energy larger 
than the ionization energy of the ion, so 7 9 is 
positive. The recombination rate is primarily de- 
termined by T\ since it depends on some inverse 
power of T (see, e.g., Arnaud and Rothenflug [3]). 
So, 7<* is negative and less than 1. It is small 
compared to 7 9 and also varies less, which fol- 
lows from evaluating Equations (6.12) and (6.15). 
(That 7 9 is larger 7 a makes physical sense, since 
7 9 is strongly affected by the high temperature 
component where 7<* is only perturbed by the 
high temperature component.) 
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and will be taken from Arnaud and Raymond 
[2] or Arnaud and Rothenflug [3] at temperature 
T\. From Equations (6.13) and (6.14), the ra- 
tio of the ionization stage densities for the bi- 
Maxwellian case can be written in terms of the 
single Maxwellian case as 


Figure 6.3. The log of the theoretical ratio of 
intensity of Fe XVI 335 A/Fe XVI 417A vs. log 
of the temperature is approximately linear over a 
broad range of temperature. The plotted curves 
are given for two sets of ionization balances, as 
discussed in Section 6 . 2 . 


(ip 1 ) • <6 - 18) 

Nx,z 1 + a 7 ar(Ti,i 2 ) V Nx,z ) Tl 

Since | 7 0 |<< 7 , ( 7 „ « -0.8, while 7 , > 10), 
Equation (6.18) can be written 


Nx^z+i , / A r x,z+i \ 
N x ,z \ Nx,z ) Tl 


1 +a 7 ,(T 1 ,T 2 ). (6.19) 


The ion densities in Equation (6.19) can be 
replaced by the spectral line intensities as in 

Ix,z+\ , ( Ix,z+i \ A' x , z +i , ( Ax.z+i \ 

Ix,z / \ Ix,z ) Tl N x ,z \ A T x,z ) Tl ’ 

( 6 . 20 ) 

since to first order the emissivities will cancel. 
(For example, for Fe XVII/Fe XVI this ratio is 
incorrect by 6 %, using calculations of Brickhouse 
et al. [7] and an a = 0.05). 

The line-ratio temperature diagnostic, for 
which a typical example is shown in Figure 6.3 
can often be fitted over a small range of temper- 
ature by 


For most cases C depends only on the two ioniza- 
tion stages involved while D depends on both the 
pair of ionization stages and the particular spec- 
tral lines used. Both C and D are determined 
from theoretical line ratios and assumed not to 
be dependent on a. The effect of the addition 
of excess high energy electrons is to increase the 
intensity ratio , so the apparent tempera- 

ture will also increase. If we subtract Equation 
(6.21) evaluated with a pure Maxwellian distri- 
bution at temperature T\ from Equation (6.21) 
evaluated with the two-component electron dis- 
tribution with line-ratio temperature T, we get 


Log 



C [Log(T) - Login )], 


( 6 . 22 ) 


since D is canceled. 

Equations (6.19) and (6.20) can be substituted 
into Equation (6.22) to produce 

Log[ 1 + Q 7 g) ^ LL [ Log[T ) — Log{T\)]. (6.23) 


Log 



= C Log(T) + D. 


( 6 . 21 ) 


Next, we use the relationship Log(x) = 
Ln(x)/Ln{ 10) and expand Ln( 1 + x) « x for 
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x « 1. For small 07,, we get 

ai q /Ln(\0) « C [Log(T) - Log(T x )], (6.24) 

which can be rewritten as 

a/3 = \Log(T) - Log(Ti ) ], (6.25) 


where 


C Ln( 10)' 


(6.26) 


Now a can be determined from two different sets 


of line-ratio temperatures (since T\ is the same 
for emission from the same plasma but the mea- 
sured temperature T will differ depending on /?) 

by _ _ 

_ LogjTj) - LogjTj) 

° A - A 


(6.27) 


where the subscripts i and j indicate the different 
sets of line ratios. 


In Table 6.1, (3 is given for different ionization 
stages using Equation (6.26) and C is evaluated 
by fitting the theoretical ratio of intensities as a 
function of the temperature. In Table 6.2 the (3 
of the two iron transitions is subtracted from the 


(3 for nickel. In both cases the values and the un- 
certainties were determined for kTi = 190±20eV 
and kT 2 = 1000±200eV. T\ is chosen to give two 
distinct components (a larger value for T 2 does 
not significantly alter the result, as will be shown 
in Figure 6.4). The line-ratio temperatures given 
in Table 6.3 are determined for both IBALl and 


IBAL2, using the intensities reported by Thomas 
and Neupert [54]. Then Equation (6.27) is used 
to determine the mixing ratio. The mixing ra- 
tio in all cases is within the uncertainty, and a is 
consistent with 0 or a small positive number in 
all cases. The assumption that Ln( 1 ± a) « a is 
invalid for large a and this would cause a to be 
underestimated by up to a factor of 2. Still, the 
observations indicate the mixing fraction is less 
than 0.25. 


Multiple temperatures along the line of sight 
could also produce a similar effect, making the 
mixing fraction lower. Since a stronger upper 
limit on the mixing fraction will be determined 


Table 6.1. Values of Parameter (3 


Lower Upper 

Stage Stage IBALl IBAL2 

Ni XVII Ni XVIII 2.58 ± 1.26 2.58 ± 1.26 

Fe XV Fe XVI 1.57 ± 0.62 1.38 ±0.66 

Fe XVI Fe XVII 1.45 ± 0.63 1.32 ±0.53 

Notes: f3 and the uncertainty in / 3 were deter- 
mined from assuming fcT\ = 190 ± 20eV, kT 2 = 
1000 ±200eV. IBALl and IBAL2 defined in Ta- 
ble 4.1. 

in Section 6.3 the effect of multiple temperatures 
along the line of sight was not investigated. The 
mixing fraction derived from the observed line 
ratios should be regarded as an upper limit on 
any bi-Maxwellian distribution. 

In Table 6.1, the j3 for Fe XV/Fe XVI is larger 
than for Fe XVI/Fe XVII. This occurs because 
C increases more than 7 9 , so the ratio actually 
decreases slightly in Equation 6.26. For a given 

а, according to Equation 6.25 the Fe XV/Fe XVI 
line-ratio temperature would be larger than the 
Fe XVI/Fe XVII line- ratio temperature, which 
is counter intuitive. The normal expectation 
would be that the measured line-ratio tempera- 
ture from the set of higher ionization levels would 
be higher. Actually the increase of the observed 
intensity ratio would be higher but the increase 
of the line-ratio temperature may be less depend- 
ing on C. 

б. 3 Theoretical Limit 

A theoretical upper limit on the mixing fraction a 
can be determined from the case where the cool- 
ing of the high energy electrons is the sole heat 
source of the corona. Then the total radiative 
flux emitted by the corona is equal to the energy 
released from the cooling of the high energy elec- 
trons. This should give a reasonable estimate of 
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Table 6.3. Temperature 


Ratio 

Fe XVI 335.400/Fe XV 417.500 
Fe XVI 360.800/Fe XV 417.500 
Fe XVII 350.500/Fe XVI 335.400 
Fe XVII 350.500/Fe XVI 360.800 
Ni XVIII 292.000/Ni XVII 249.200 


IBAL1 IBAL2 
6.45 ±0.09 6.41 ±0.07 

6.43 ± 0.09 6.39 ± 0.07 
6.50 ±0.07 6.37 ±0.05 
6.52 ± 0.08 6.38 ± 0.06 
6.45 ±0.14 6.45 ±0.14 


Notes: The line-ratio temperatures were calculated using the intensities given by Thomas and 
Neupert [54], assuming Landini and Monsignori-Fossi [41] emissivities. The uncertainties includes 
the fitting uncertainty, and the instrument uncertainty described by Thomas and Neupert. The 
uncertainty also includes the ionization balance uncertainty and the emissivity uncertainty. 


Table 6.4. Mixing Fraction, a 


Ratio IBAL1 

Fe XVI 335.400/Fe XV 417.500 0.00 ± 0.16 

Fe XVI 360.800/Fe XV 417.500 0.02 ± 0.16 

Fe XVII 350.500/Fe XVI 335.400 -0.04 ± 0.14 

Fe XVII 350.500/Fe XVI 360.800 -0.06 ± 0.15 


Table 6.2. Values of Parameter (/3 Ni - fi Fe ) 


Lower Upper /?/v, — (3 Fe 

Stages Stages IBAL1 IBAL2 
Fe XV Fe XVI 1.02 ±0.64 1.20 ±0.71 

Fe XVI Fe XVII 1.13 ± 0.64 1.26 ±0.68 

Notes: The uncertainties in the difference in (3 
were determined assuming kTi = 190 ± 20eV, 
kT 2 = 1000 ± 200eV. (The uncertainty is less 
than would be determined by taking the square 
root of the sum of the squared uncertainties in 
Table 6.1 since they are not independent.) 


IBAL2 
0.03 ±0.13 
0.05 ±0.14 
0.06 ±0.12 
0.06 ±0.12 


the mixing fraction if conduction, which is ne- 
glected, is not too important and if the plasma 
is in a steady state. 

The time for equipartition for two sets of 
Maxwellian particles with different temperatures 
is given by Spitzer (1961) [51] as 


A\ ' A 2 J ’ 

(6.28) 

where A is the ratio of the Debye shielding dis- 
tance to the impact parameter, T, is the temper- 
ature of distribution i, A{ is the atomic mass of 
the particles of group i and N\ is the bulk plasma 
density. For the case where both sets of particles 
are electrons and 

N t =N e (l -a), (6.29) 

we get 

‘^^N^T) lT ' + T ^ ec - (6 ' 30) 


A \ iT-2 

f K 07 1 1 

teq - Nl ZlZ*ln{A) 
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(For electrons with N x = 10 9 cm -3 , T x = 2 x 
10 6 /\, and T 2 = lO 7 /^ we get « 0.3sec.) We 
take T 2 and N 2 to be the temperature and density 
of the high-temperature component and use 


which can be rearranged as 


^ 5.87eA°- 5 (T x + r 2 ) 15 
ln(\)k(T 2 — T x ) ’ 


(6.39) 


N 2 = N e a. 

Then the exponential decay rate is 


(6.31) 


t = („(!- a) (6.32) 


if T x and N x T x -\-N 2 T 2 are constant (steady state). 

Each electron which is transferred from 
the high-temperature component to the low- 
temperature component would release an amount 
of energy given by k(T 2 — T x ). Thus, the total 
heating of the bulk plasma in steady state would 
be 


F h = 


k{T 2 — T\)N 2 


(6.33) 


Substituting Equations (6.30) and (6.32) in 
Equation (6.33) we get 


For coronal conditions /n( A) is nearly con- 
stant. According to tables in Spitzer (1961) on 
pg 128 [51], /n(A) = 19.3 for T = 10 6 A' and N e = 
10 9 /cm 3 . Cook [9] gives £ = 2.35 x 10 ~ 23 erg cm 3 
for log{T ) = 6.4, and e = 6.735 x 10 ~ 23 erg cm 3 
for log(T) = 6.3, using Meyer 1985 [37] elemen- 
tal abundances. (If Meyer 1992 [40] composition 
is used, then the emissivities, and consequently 
the upper limits on a, will increase by roughly 
a factor of 3.5.) Thus, using Meyer 1985, for 
log(T\ ) = 6.4, we get 


a - 


i.2i x io- 9 (r! + r 2 ) 15 

t 2 -t x 


and, for log(T x ) = 6.3, 


(6.40) 


k{T 2 - T x )N x N 2 ln( A) 
5.87A°- 5 (T 1 +T 2 ) 15 (1 - a)' 


(6.34) 


a = 


3.46 x 10- 9 (7\ + T 2 ) 15 

t 2 -t x 


(6.41) 


Replacing N x and N 2 in Equation (6.34) then 
gives 


k{T 2 - T x )N 2 e a ln{ A) 
5.87 A°- 5 (T X + T 2 ) 15 ' 


(6.35) 


The radiative flux, discussed in Section 5.2, is 
given by 

F r = eN 2 . (6.36) 

Where the plasma emissivity is primarily depen- 
dent on the bulk plasma temperature T x . If equi- 
librium exists and conduction is not important, 
then heating for the high-energy electrons needs 
to be equal to or less than the radiative flux. So 
an upper limit on the mixing fraction, a, can be 
determined, by setting F h > F r yielding 


or 


eN i > HT, - 7 \) Nla ln( A) 


£ > 


5.87 A°- 5 (T X + T 2 ) 15 ’ 

k(T 2 -T x )aln( A) 
5.87A°- 5 (T 1 +T 2 ) 1 - 5 ' 


(6.37) 

(6.38) 


Both curves for a are plotted in Figure 6.4 versus 
T 2 . As can be seen in the figure, for T 2 > 10 7 K 
the upper limit on the mixing fraction is rela- 
tively flat. Any additional heating mechanism 
would cause a to be smaller. If heat conduction 
is included the upper limit could be increased by 
about a factor of 2, but it will still be small for 
the steady state case 

The theoretical upper limit on a of less than 
10 -5 does indicate that the previous results 
Maxwellian assumption is reasonable for Chap- 
ters 2, 4 and 5. With a typical 0 of 20 the temper- 
ature changes by less than one tenth of a percent, 
which is small compared to the other errors. 

Since the upper limit on a is small, trying to 
observe it will be extremely difficult experimen- 
tally. The best possibility will be to look at the 
ratio of intensities between helium-like ionization 
stages and hydrogen-like ionization stages where 
the ionization energy is high. For example Si 
XIII is helium like and has an ionization energy 
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Figure 6.4. The upper limits on a as a function of 
T 2 for log(T\ ) = 6.3 and log{T\) — 6.4 are shown. 
This assumes that heating done by the high- 
temperature component is radiated away with no 
other heat source, and that conduction is negli- 
gible 


of 523.52 eV, while Si XIV has an ionization en- 
ergy of 2437.76 eV according to Cowan [10]. Ac- 
cording to Arnaud and Rothenfiug [3] 10 -0 ' 24 of 
the silicon should be Si XIII while 10 -5 97 should 
be Si XIV. Both have spectral lines between 4- 
7A. Using kT\ = 190eV and kT 2 — 1000eV the 
amount of Si XIV will double for an a — 10~ 3 . 
This effect will increase for elements with higher 
Z elements. In the case of S XV and S XVI the 
effect will be similar for a = 10 -4 . For lower Z 
elements such as Ne IX to Ne X, an a — 0.1 will 
approximately double the amount of Ne X. The 
effect of a is stronger on the higher Z elements, 
at any given temperature 7\. The trouble is that 
often the hydrogen-like lines are very weak com- 
pared to the helium-like lines and will be hard to 
detect. Any significant quantity of 10 7 A' plasma 
along the line of sight would make the analysis 
impossible, since it will be the stronger emitter 
of the Si XIV and S XVI spectral lines. 
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Chapter 7 
Conclusions 


7.1 Relative Elemental 
Abundances 

The present understanding of coronal elemental 
abundance is discussed in Section 2.2. In Sec- 
tion 4.6 the set of relative elemental abundances 
for silicon, aluminum, chromium, and iron was 
determined. The relative elemental abundance 
for silicon to iron for SERTS 1989 using IBALl 
agrees within uncertainties with either photo- 
spheric (Allen [1]) or coronal elemental abun- 
dances given by Meyer [38, 40]. 

However, The relative elemental abundance for 
SERTS 1991 differs from elemental abundances 
given in the literature for both photospheric and 
coronal elemental abundances. For SERTS 1991, 
the lower the First Ionization Potential (FIP) 
of the element the larger the enhancement of 
the element compared to photospheric or coronal 
abundances. A tilt in the upper level of the step 
function in Figure 2.1 would explain this result. 

The differences between SERTS 1989 and 
SERTS 1991 sets of elemental abundances indi- 
cate that the solar coronal elemental abundance 
can vary. A simple step function would not ex- 
plain all solar coronal conditions. Both Saba and 
Strong [48] and McKenzie and Feldman [34] re- 
ported variations in the relative elemental abun- 
dances for low FIP elements (Section 2.2). Ad- 
ditional observations might determine how best 
to model the coronal elemental abundance. 

The difference in the SERTS 1989 and 1991 el- 


emental abundance is real. The Si XI 303 A, Fe 
XVI 361 A, and the Fe XIV 334 A spectral lines 
were used for both flights. The measured line- 
ratio temperatures were similar. Any change in 
the ionization balance calculations would change 
the Fe to Si abundance for both flights by simi- 
lar amounts. The difference in relative elemental 
abundance between the two flights would still re- 
main. 

The most likely explanation would be that the 
observed regions of the solar corona were dif- 
ferent. Both were observed off the solar limb. 
Both were quiet solar regions with some active 
regions near the line of sight. The primary differ- 
ence between the two flights was that for SERTS 
1991 the emission measure was a factor of ten 
times larger than for SERTS 1989, a hot loop 
was nearby, and the Sun was more active near the 
line of sight. The two flights were 2 years apart so 
temporal variations of coronal abundances could 
also explain the result. 

7.2 Radial Dependence of 
Coronal Parameters 

The three observation sites examined in this dis- 
sertation all have line-ratio temperatures which 
increase with height, where as the line-ratio den- 
sities, emission measures, and spectral line inten- 
sities decrease with height above 1.03 R®. The 
decrease of the emission measure is roughly expo- 
nential. The line-ratio temperature increase ap- 
pears to be linear, with any second-order depen- 
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dence being smaller than the uncertainty in the 
measurement of its dependence. These smoothly 
varying radial dependencies allowed the use of 
simple coronal models with spherical symmetry. 

The line-ratio temperatures and scale-height 
temperatures agree to within approximately 25% 
(Table 4.6). The scale-height temperature was 
derived from the radial dependence of the inten- 
sity. The line-ratio temperature was determined 
from the ratio of the intensities of different spec- 
tral lines. These two temperature measurements 
were nearly independent. Since the line-ratio 
temperature requires fewer assumptions and can 
be used to determine the radial temperature 
dependence, it is the preferred method for de- 
termining the coronal temperature. Reliable 
determination of the line-ratio temperature re- 
quires many spectral lines from different ioniza- 
tion stages of the same element, where none of 
the lines is strongly affected by resonance scat- 
tering. 


7.4 Bi-Maxwellian 

For a bi-Maxwellian electron distribution to af- 
fect the ionization balance, we have shown (Sec- 
tion 6.3) that, under steady-state conditions, the 
mixing fraction must greatly exceed the value of 
2 x 10~ 5 that we determined to be an upper limit 
for the corona, assuming that the cooling of high 
energy electrons is the sole source of energy ra- 
diated away. Therefore, any bi-Maxwellian mix- 
ing fraction that can exist under solar coronal 
steady-state conditions will not affect the ioniza- 
tion balance. 


7.3 Heating 

In all cases, the line-ratio temperature increases 
with height. To supply the heat which would 
be conducted downward due to the tempera- 
ture gradient and radiation, a source of heat is 
required at some location above the maximum 
observed height (1.15 R@ for SERTS 1989 and 
1.2 R q for SERTS 1991). 

The maximum divergence of heat flux was de- 
termined from the radial dependence of the tem- 
perature. The total radiative power was obtained 
from the emission measure and line-ratio temper- 
ature (Section 5.2). By comparing the divergence 
of heat flux and the total radiative power, for a 
filling factor of 1.0, heating was shown to be nec- 
essary for two of the three cases over some of the 
region below 1.15 /?©, with the third case being 
inconclusive. For a filling factor of 0.1, heating is 
necessary for all three cases over the entire region 
observed. 
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Chapter 8 


Future Applications of SERTS and 
SOHO 


8.1 SERTS 

Additional off-the-limb SERTS observations are 
highly desirable for the study of the coronal el- 
emental abundances and heating requirements. 
Use of a narrower lobe (because there would be 
fewer overlapping emission lines) should increase 
the number of usable lines on both sides of the 
lobe. Additional iron spectral lines from different 
stages of ionization should lead to a better de- 
termination of the temperature or temperature 
profile along the line of sight. This would al- 
low stronger constraints to be placed on heating. 
Additional spectral lines from different elements 
would allow better determination of the relative 
elemental abundances. 

Joint observations from Yohkoh and other 
space platforms would allow a more extensive 
analysis of coronal temperature structure and 
composition. The Yohkoh satellite observes hot- 
ter plasma along the line of sight. Coronagraph 
measurements would allow large-scale coronal 
structures to be studied. 

8.2 SOHO CDS/NIS 

The type of analysis done in this dissertation 
could also be applied to data obtained from the 
planned SOHO satellite (Solar and Heliospheric 
Observatory). For example, using the 4 by 240 
arcsec slit of the SOHO CDS/NIS instrument 


(Coronal Diagnostic Spectrometer/Normal Inci- 
dence Spectrometer), and rastering the slit off 
the limb to form images from a number of spec- 
tral lines, would provide a useful data set for this 
kind of analysis. The 4 by 240 slit would have 
a resolution of 0.168 A, which is larger than the 
SERTS narrow slit, and is smaller than the lobes. 
Therefore, more lines would be usable than were 
used in this study. This would allow better de- 
terminations of emission measure, density, and 
temperature. These results would allow the de- 
termination of the relative elemental abundances 
and imposing of further heating constraints. 

The addition of more spectral lines would al- 
low a differential emission measure to be deter- 
mined. A differential emission measure deter- 
mines the emission measure distribution over a 
range of temperatures by using emissivity of the 
different spectral lines as a function of tempera- 
ture, and finding the best overall fit to the data. 
This was not done for this dissertation, because 
there were insufficient resolvable lines from the 
same element to give a well constrained result. 
For the cases studied, an uniform temperature 
gave a reasonable fit to the data. 
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Appendix A 

Solar Plasma Parameters 


Region 

Photosphere 

Chromosphere 

Corona 


Temperature 
6,000K 
10 A K 
2 x 10 6 K 


Density 

H s 

u p 

10 16 cm~ 3 

300km 

2n x 10 12 Hz 

10 n cm -3 

500km 

67r x 10 9 Hz 

10 9 cm -3 

90,000km 

6tt x 10 s Hz 


Useful Equations 


Plasma Frequency, ~ 10 4 N 05 (rad/ sec) 

Mean-Free-Path for fully ionized plasma, A m /p « 3.8 x 10 9 ^i m 
Hydrostatic Scale Height, H„ w i 
Temperature 10 4 A" ~ 1 eV 

Time for Equipartition for Two Sets of Maxwellian Particles with Different Temperatures (Equa- 
tion 6.28) 


t 


eq 


5.87 


^ 1^2 

7V 1 Z, 2 Z 2 2 /n(A) 



(A.l) 


Example Ni — 10 9 cm~ 3 , T\ = 2 x 10 6 A , and Ti — 10 7 A' 
For Electrons (t eq & 0.3 sec) 

For Protons ( t eq « lOsec) 

For Protons (Ti) and Electrons (T 2 ) ( t eq « 10 minutes) 
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Glossary of Symbols 


Symbol 

Description 

Page First 



used in 

A 

Area 

9 

A, 

Mass of particle i in amu 

76 

a 

2nd fitting parameter for line-ratio temperature 

41 

C 

Correction factor from line intensity equation 

17 

F(E,T ,) 

Maxwellian Electron Distribution 

72 

F c 

Electron heat conduction 

66 

F' c 

Divergence of the heat conduction 

66 

F h 

total heating Flux 

77 

F 0 

Constant part of the divergence of the heat conduction 

66 

F T 

Radiative power 

66 

Ft 

Electron Distribution 

72 

G 

Universal constant of gravity (6.67 x 10 ~ s dyn cm 2 g~ 2 ) 

9 

9 

Acceleration due to gravity 

9 

Has 

Observed density scale height 

21 

H, 

Density scale height 

10 

R's 

Effective density scale height 

10 

k 

Boltzmann constant(8.617 x 10 _5 eF A' -1 ) 

10 

M e 

Solar mass (1.989 x 10 33 gm) 

9 

N e 

Electron density 

10 

Ne 

Average of the electron density squared 

17 

Ni 

Density of Particle Distribution i 

76 

p 

Pressure 

9 

Rq 

Solar Radius (6.9599 x lO 10 cm) 

10 

r 

Radius (distance from Sun Center) 

9 

To 

Reference radius (Section 2.1.1) 

10 

r’ 

Line-of-sight radius 

17 

T 

Temperature 

10 

T e 

Electron temperature 

13 

T t 

Temperature of Particle Distribution i 

76 

T 

1st fitting parameter for line-ratio temperature 

20 

Lg 

time for equipartition 

76 

y 

(i-p) 

21 

z 

Length along the magnetic field line 

66 

Zt 

number of protons of Particle i 

76 
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Symbol 

Description 

Page First 



used 

in 

7 q 

Correction Term for Ionization Rate 

73 


la 

Correction Term for Recombination Rate 

73 


£T 

Total plasma emissivity 

66 


c 

Correction factor from line intensity equation 

17 


A 

the ratio of the Debye shielding distance 

76 



to the impact parameter 



^mfp 

Mean Free Path 

71 


P 

Average atomic weight 

10 


P 

Mass density of the gas 

9 


0 

Conversion Factor 



Ix,z 

Spectral line intensity 




(Following symbols are defined along line of 

sight) 


Symbol 

Description 


Page First 




used in 

/ 

Filling factor 


17 

L 

Effective line-of-sight length 


17 

N e 

Line-of-sight density 


16 

Qx,z,i,j 

Density weighted average contribution function 

16 

T e 

Line-of-sight temperature 


16 

X 

Path length along the line of sight 


16 

1 

Column emission measure 


17 
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For the following, X=element, Z= ionization stage, 
i is an excitation level, and j is a second excitation level 


(Following symbols are defined along line of sight) 


Symbol 

Description 

Page First 
used in 

A x,z,i,j 

Einstein spontaneous emission rate 

13 

b x,z,i,j 

Branching ratio 

15 

Cij 

Electron collision excitation rate (j > i ) 

14 

Cj4 

Collision de-excitation rate (j > i ) 

14 

Ex,z,i 

Energy of an excitation level 

14 

E x,z,i,j 

Energy per photon 

13 

Ix,Z,i,j 

Spectral line intensity 

15 

N x,z,i 

Density 

13 

qx,z 

Ionization rate coefficient 

15 

Qx,Z,i,j 

Contribution function for a spectral line 

13 

a x,z 

Radiative recombination rate coefficient 

15 

Iji 

Effective collision strength 

14 

e x,z,i,j 

Emission per unit volume 

13 


Statistical weight for level j 

14 


NOTES: 

dummy indices not listed 
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mu in i 


The general notation used to indicate the spectral line or atomic energy level in question needs to 
be introduced. Since this thesis compares lines from different elements, ionization stages, levels of 
excitation, and solar radial positions, the notation needs up to four subscripts. An example of the 
notation is given by EFe,xv,s, 1 , where E indicates the energy per photon, Fe is the element, XV is 
the ionization stage, 5 denotes the initial level and 1 the ground state too which the ion de-excites. 
The general case of an arbitrary emission line produced by an arbitrary transition in an arbitrary 
element would be written Ex,z,i,j where X is the element and Z is the ionization stage. Densities 
are specified with three subscripts or less, for example with NF e ,xv ,5 indicating the volume density 
of Fe XV in the fifth excitation level. The density of an element in a particular ionization stage, 
summed over all excitation levels will be Nx,Zi where 

Nx,z = Nx,z,j- (A-2) 

j 

Finally radial position is denoted in the common way as an argument in parenthesis. For example 
Nx{r) is the volume density of element X at radius r. 

(Following symbols are defined along line of sight) 

Symbol Description Page First 

used in 


Different 

Qx,z Ionization Cross Section 73 

E Kinetic Energy of the Particle 71 

a Mixing Fraction 72 

C Fitting Parameters 72 

D Fitting Parameters 72 

(Tj) Measured Line-ratio temperature of line-ratio i 74 
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